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When a jet collides coherently with many parton scatterers at very high energies, the Bose-Einstein 
symmetry with respect to the interchange of the virtual bosons leads to a destructive interference of 
the Feynman amplitudes in most regions of the momentum transfer phase space but a constructive 
interference in some other regions of the momentum transfer phase space. As a consequence, the 
recoiling scatterers have a tendency to come out collectively along the incident jet direction, each 
carrying a substantial fraction of the incident jet longitudinal momentum. The manifestation of 
the Bose-Einstein interference as collective recoils of the scatterers along the jet direction may have 
been observed in angular correlations of hadrons associated with a high-pT trigger in high-energy 
heavy-nuclei collisions at RHIC and LHC. 
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I. INTRODUCTION 

Recently at RHIC and LHC, angular correlations of 
produced hadron pairs in AuAu, PbPb, and pp colli- 
sions have been measured to obtain the yield of produced 
pairs as a function of A(j) and A77, where A0 and A77 are 
the azimuthal angle and pseudorapidity differences of the 
produced pair. The correlations appear in the form of a 
"ridge" that is narrow in A<f> at A(/) ~ and Acf) ~ tt, but 
relatively flat in A77. They have been observed in high- 
energy AuAu collisions at RHIC by the STAR Collabora- 
tion [M^, the PHENIX Collaboration H^-Hl, and the 
PHOBOS Collaboration [25|, with or without a high-pT 
trigger 0, [l^ - fisj . They have also been observed in pp 
and PbPb colhsions at the LHC by the CMS Collabo- 
ration [H-ii], the ATLAS CoUaboration and the 
ALICE Collaboration 

Subsequent to the observation of the A6-Ar] angular 
correlations, a momentum kick model |3ll439f was put 
forth to explain the "r idg e" phenomenon, along with 
many other models |40l - [74| . The model identifies ridge 
particles as medium partons, because of the centrality 
dependence of the ridge particle yield and the similar- 
ity between the temperature and baryon/mcson ratio of 
ridge particles with those of the bulk medium. The model 
assumes that these medium partons have an initial ra- 
pidity plateau distribution, and when they suffer a col- 
lision with a jet produced in the collision, they receive 
a momentum kick along the jet direction, which will be 
designated as the longitudinal direction. The physical 
contents of the momentum kick model lead to features 
consistent with experimental observations: (i) the lon- 
gitudinal momentum kicks from the jet to the medium 
partons along the jet directiongive rise to the A<f) ~ 
correlations on the near-side [l|,ll[3ll-l33, (ii) the medium 
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parton initial rapidity plateau distribution shows up as a 
ridge along the Ar] direction P, S UMl 113, H!] , ( iii) the 
displaced parton momentum distribution of the kicked 
medium partons leads to a peak at pr^ 1 GeV/c in the 
Pt distribution of the kicked partons [S, [SI, UMl : and (iv) 
the kicked medium partons therefore possess a correla- 
tion of IpitI ~ |p2tH1 GeV, without a hi gh-p r trigger, 
for both the near side and the away side [iSl, l38[ . The 
model has been successful in describing extensive sets of 
triggered associated particle data of the STAR Collab- 
oration, the PHENIX Collaboration, and the PHOBOS 
Collaboration, over large regions of pt, Acj), and Arj, in 
many different phase space cuts and pt combinations, 
including dependencies on centralities, dependencies on 
nucleus sizes, and dependencies on collision energies 

The phenomenological success of the momentum kick 
model raises relevant questions about the theoretical 
foundations for its basic assumptions. While the rapidity 
plateau structure of the medium parton distribution may 
have its origin in the Wigner function of particles pro- 
duced in the fragmentation of a string (or flux tube) [3al , 
or in a color-glass condensate [l^l, the origin of the pos- 
tulated longitudinal momentum kick along the jet direc- 
tion in the model poses an interesting puzzle. If each of 
the jet-(medium parton) collision were a two-body elastic 
collision, it would lead only to a dominantly forward scat- 
tering and a small longitudinal momentum transfer to the 
medium scatterers along the jet direction. The hard lon- 
gitudinal momentum kick (of order GeV) received by the 
medium scatterers along the jet direction in the momen- 
tum kick model is unlikely to originate from incoherent 
two-body collisions. 

From an observational viewpoint, the experimental ob- 
servation that many medium partons come out along the 
jet longitudinal direction with A<^ ~ implies the pres- 
ence of collective recoils of the medium partons, each 
of which must have acquired a substantial longitudinal 
momentum transfer (momentum kick) along the jet di- 
rection from the jet. What is the origin of such collective 
recoils, in the interaction of the jet with medium partons? 
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As the high-pT trigger particles under experimental 
consideration have a momentum of order a few GeV 
H, [20I423I [27l - [30j , we therefore consider incident jets with 
an initial momentum of order 10 GeV. Strictly speaking, 
jets with an initial momentum of order 10 GeV should be 
more appropriately called "mini-jets" [t^. For brevity 
of nomenclature, we continue to use the term "jet" to 
represent mini-jets. 

In the collision of a jet p={po,p) with n medium par- 
tons, how dense must the medium be for the multiple 
collisions to become a single coherent (l-l-n)-body colli- 
sion, instead of a sequence of n incoherent 2-body colli- 
sions? It is instructive to find out the conditions on the 
medium density and the energy of the incident jet that 
determine whether the set of multiple collisions are co- 
herent or incoherent. For such a purpose, we consider 
a binary collision between the incident fast particle p 
and a medium scatterer with the exchange of a bo- 
son, in the medium center-of-niomentum frame. The 
longitudinal momentum transfer Qz for the binary col- 
lision can be obtained from the transverse momentum 
transfer qt by Qz ^ q^/2pQ. The longitudinal momen- 
tum transfer is associated with a longitudinal coherence 
length Azcoh ^ ^/'Zz ^ '^^Po/Qt that specifies the uncer- 
tainties in the longitudinal locations at which the virtual 
boson is exchanged between the fast particle and the scat- 
terer. The time it takes for the fast particle to travel the 
distance of the longitudinal coherence length, 2poh./q^c, 
can also be called the virtual boson formation time. 

The nature of the multiple scattering process can be 
inferred by comparing the longitudinal coherence length 
Azcoh with the mean free path A of the jet in the dense 
medium that depends not only on the density of the 
medium but also on the binary collision cross section. 

If Azcoh ^ A, then a single binary collision is well 
completed before another binary collision begins, and the 
multiple collision process consists of a sequence of n in- 
coherent two-body collisions. On the other extreme, if 

^^coh 

- 2hpa/ql > A, (1) 

then a single binary collision is not completed before an- 
other one begins, and the multiple collision process con- 
sists of a set of coherent collisions as a single (l-l-n)-body 
collision. For a set of initial and final states in such a co- 
herent (l-l-n)-body collision, there are n\ different trajec- 
tories in the sequence of collisions along Azcoh at which 
various virtual bosons are exchanged. By Bose-Einstein 
symmetry, the total Feynman amplitude is then the sum 
of the n! amplitudes for all possible interchanges of the 
exchanged virtual bosons. The summation of these n\ 
Feynman amplitudes and the accompanying interference 
constitute the Bose-Einstein interference in the passage 
of the fast particle in the dense medium. 

In high energy central collisions between heavy nuclei 
such as those at RHIC and LHC, both jets and dense 
medium are produced after each collision. The jets will 
collide with partons in the dense medium, and these col- 
lisions may satisfy the condition for coherent collisions. 



In a binary collision at RHIC and LHC, the longitudi- 
nal coherent length Azcoh is of order 25 fm, for a typical 
transverse momentum transfer of gT'~-"0.4 GeV/c from a 
jet of momentum po~10 GeV/c to a medium parton. The 
longitudinal coherent length Azcoh is much greater than 
the radius i? of a large nucleus. On the other hand, the 
away side jet is quenched by the dense medium in the 
most central AuAu collisions at RHIC and LHC [zl-il, 
and the near-side jet collides with about 4-6 medium par- 
tons [3l[-[3^. Therefore, one can infer that the mean- free 
path A for the collision of the jet with medium partons is 
much smaller than the nuclear radius R. In high-energy 
central nuclear collisions at RHIC and LHC, the collision 
of a fast jet with medium scatterers satisfy the following 
condition 

Azcoh - 2hpo/qT > i? > A. (2) 

As a consequence, the multiple collision process consti- 
tutes a set of coherent collisions. There will be Bose- 
Einstein interference effects in the passage of the jet in 
the dense medium. 

In this connection of Bose-Einstein interference, we 
note that Bose-Einstein interference effects in high- 
energy QED and QCD collisions have been observed pre- 
viously by many workers [83l - l90j . In the emission or ab- 
sorption of n identical bosons from an energetic fermion 
in Abelian and non-Abelian gauge theories leading to an 
on-mass-shell final fermion, the Bose-Einstein symmetry 
with respect to the interchange of virtual bosons leads to 
the sum of a the set of n\ Feynman amplitudes that turns 
out simply to be a product of delta functions. These 
distributions may be thought of as peaked interference 
patterns produced by the coherent addition of various 
symmetrized broad Feynman amplitudes [o^l ■ In the col- 
lision of two fermions, the sum of the ladder and cross- 
ladder diagrams also exhibits remarkable Bose-Einstein 
interference leading to similar products of delta functions 
and the eikonal approximation [sst - fooj . 

We would like to study similar interference effects for 
the case of the multiple collisions of a jet parton with 
many parton scatterers. We wish to explore whether 
such a Bose-Einstein interference may be the origin of the 
longitudinal momentum kick. As the interference arises 
from the coherent collision of more than one particle, it is 
clearly a many-body effect that occurs with two or more 
scatterers. 

The interference of Feynman amplitudes is only one 
of the effects of coherent collisions. There is another 
important effect that accompanies the coherent 
body collision and changes the nature of the collision 
process. It shows up as an increase in the number of 
degrees of freedom. In a sequence of n incoherent two- 
body collisions or quasi-elastic collisions, there are two 
degrees of freedom in each two-body collision, which can 
be chosen to be q^rp = (qix,(j>i). The sequence of n 
incoherent or quasi-elastic two-body collisions contains 
only 2n transverse momentum transfer degrees of free- 
dom, {qr^j,, qr2Tj 93Ti 9raT}- The longitudinal momen- 
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turn transfer g.j^ in each individual jet-parton collision 
is a dependent variable, depending on the correspond- 
ing transverse momentum transfer as qiz |qjyp/2po- 
In contrast, in the case of coherent collisions of the jet 
with n particles, there are 3n — 1 degrees of freedom, af- 
ter the 3(n -I- 1) degrees of freedom are reduced by the 
constraints of the conservation of energy and momen- 
tum. The degrees of freedom increases from 2n for inco- 
herent collision and quasi-elastic collisions to 3n — 1 for 
coherent collisions. We can choose the 3n — 1 indepen- 
dent variables to be {q'lT'Ziz, 92t92z, Qst^Sz, 9„t9"2} 
for coherent collisions, subject to a single condition of an 
overall energy conservation. Thus in the case of coher- 
ent collisions, the set of longitudinal momentum trans- 
fers, {qizji — l:2,...,n}, can also be independent vari- 
ables with their own probability distribution functions, 
depending on the sum of nl symmetrized Feynman am- 
plitudes. The recoiling scattcrers share the longitudinal 
momentum of the incident jet. The Qiz degrees of free- 
dom allow the medium partons to acquire substantial 
fractions of the longitudinal momentum of the incident 
jet, as we shall demonstrate in Section V of this paper. 

In the problem of the passage of a jet in a dense 
medium, the multiple collision process of the jet with 
medium partons is usually examined in the potential 
model [Za-Ill as in the Glauber theory [9|. It 

should be noted that the occurrence of coherent collision 
discussed here coincide with the condition for multiple 
scattering in the Glauber theory given by Eq. (173) of 
[al as 

Poa^/h > R, (3) 

where R is the dimension of the medium, and a is the 
range of the interaction which can be related to the av- 
erage transverse momentum transfer by a ~ Ti/^t- Thus 
the Galuber's condition for multiple scattering is 

fipjqr > (4) 

As the dimension of the medium is presumably much 
greater than the mean free path A between collisions, 
condition ([3]) for Glauber's multiple scattering is consis- 
tent with 

hpz/qr > i? > A (5) 

that is nearly the same as the coherent condition ([2]) we 
discussed here, except for an unimportant difFercncc of a 
factor of two. 

The near-coincidence of the conditions for coherent 
collisions of Eqs. ([2]) and (O indicates that the total- 
ity of coherent collisions can be further divided accord- 
ing to the degree of incident jet longitudinal momentum 
losses. There are elastic and quasi-elastic coherent col- 
lisions with nearly no scatterer recoils and little loss of 
the incident jet longitudinal momentum that are at the 
center of Glauber-type potential model theories. There 
are also coherent collisions in which scatterers can un- 
dergo longitudinal recoils and the incident jet can lose a 



substantial longitudinal momentum arising from the re- 
coils of scatterers. The elastic and quasi-elastic coherent 
collision processes can be adequately studied in the po- 
tential model [ttI - IsiI [9ll - l93 |. with the assumption that 
the scatterers suffer approximately no recoils and appear 
in the form of potential centers around which the projec- 
tile traverses. 

We wish to investigate coherent collisions in the entire 
domain of longitudinal scatterers recoils and the corre- 
sponding incident jet longitudinal momentum loss due 
to these recoils. For a general treatment of longitudinal 
recoils of scatterers, the Glauber-type potential model 
cannot really be used. It is inadequate because it is re- 
stricted to elastic and quasi-elastic collisions and does not 
allow the scatterer longitudinal recoils to be independent 
degrees of freedom. For well-founded reasons which we 
shall discuss in detail in Section II, we shall forgo the 
potential model but turn to the Feynman amplitude ap- 
proach to study the longitudinal recoils of the scatterers 
in a general treatment of coherent collisions. 

In the remaining parts of the manuscript, we shall con- 
sider implicitly only coherent multiple collisions of the jet 
with medium partons unless indicated otherwise. In Sec- 
tion III, we examine (coherent) multiple collisions of an 
energetic fermion with fermion scatterers in the Abelian 
gauge theory to discuss how the Bose-Einstein symme- 
try in the Feynman amplitudes can lead to a destruc- 
tive interference over most regions of phase space but a 
constructive interference in some other regions of phase 
space. In Section IV, we study the cross section for such 
a multiple collision process. In Section V, we examine 
the consequences of the constraints on the recoils of the 
scatterers and find that the scattcrers tend to come out 
collectively along the incident jet direction, each scatterer 
acquiring approximately l/2n fraction of the longitudinal 
momentum of the incident jet. In Section VI, we examine 
the case for the coherent collision of an energetic fermion 
with fermion scattcrers in the non- Abelian gauge theory. 
In Section VII, we generalize our considerations to the 
collision of a gluon jet on quark scatterers. In Section 
VIII, we examine the collision of a gluon or a fermion jet 
on gluon scatterers. In Section IX, we study the longi- 
tudinal momentum distribution for gluon scatterers and 
distinguish it from the longitudinal momentum distribu- 
tion for quark scatterers. In Section X, we examine the 
signatures of the Bose-Einstein interference in the pas- 
sage of a jet in a medium and compare these signatures 
with experimental data. In Section XI, we present our 
conclusions and discussions. 



II. THE POTENTIAL MODEL VERSUS THE 
FEYNMAN AMPLITUDE APPROACH 

The potential model of multiple collisions [ttI - IsTI . IqTI - 
Iq^ I has been a useful tool to study the energy loss of a 
jet as it passes through a dense medium in elastic and 
quasi-elastic scatterings. It can also be adequately used 
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to study the longitudinal momentum transfers qi^ from 
the jet to the scatterers in these elastic and quasi-clastic 
processes. The potential model gives the longitudinal 
recoil momentum 

- |g,rlV2po, (6) 

for each jet-scatterer collision, in either coherent or inco- 
herent collisions. 

We shall sec that in the Feynman amplitude approach 
the probability distribution for the longitudinal recoil g,^ 
includes not only the probability for quasi-elastic pro- 
cesses with little scatterers longitudinal recoils but also 
a significant probability for the scatterers to recoil with 
a substantial fraction of the incident jet longitudinal mo- 
mentum. It is therefore important to inquire which of 
the two different approaches should be used to investi- 
gate the longitudinal recoils of the scatterers in a general 
treatment of coherent collisions that are at the focus of 
our attention. 

It is instructive to review how the longitudinal recoils 
of the scatterers are determined in the two different ap- 
proaches. In the potential model [ttI - IsiI . loTI - fo^ . the 
medium parton scatterers are assumed to be fixed in 
space and the scatterings between the incident particle 
and the medium scatterers are represented by potentials 
with the positions of the scatterers as static centers, as 
in the Glauber model for potential scattering [o^l- By 
using static scatterer centers, the range of applicability 
of the potential model is limited to elastic or quasi-elastic 
scatterings with small momentum transfer <7j. Further- 
more, we can trace how qi={qiTi Qiz) is determined in the 
potential model, in the collision of the jet with the zth 
scatterer. After passing through the potential generated 
by the ith scatterer, the trajectory of the incident jet is 
obtained and the transverse momentum transfer q^rp de- 
termined. The longitudinal momentum transfer qi^ from 
the jet to the scatterer is subsequently calculated to en- 
sure longitudinal momentum conservation. Hence the 
longitudinal momentum transfer from the jet to the ith 
scatterer is determined as qiz ~ qfj^/2pQ, as given in Eq. 
([5]). As the longitudinal momentum transfer qiz is de- 
termined completely by and po, it is no longer an 
independent dynamical variable. It is in effect a frozen 
and dependent quantity. Not only is the potential model 
limited to scattering with small momentum transfers q^, 
the model is also limited to cases in which the longitudi- 
nal momentum transfers qiz from the jet to the scatterers 
are no longer independent dynamical variables. 

In the potential model with many scatterers, the total 
potential is the sum of all scatterer potentials. In the 
coherent collisions limit, the wave function of the incident 
jet becomes 

ip{z) = exp{ifc(z - zo) -I- ixtot{b, z)}, (7) 
where the phase xtot{b,z) is the sum of the scattering 



phases from scatterer potentials Vi centered at {bi,Zi}, 

Xtot{b,z)^-— dz'Y,V^ib-b^,z' - Z^). (8) 

The width of the transverse momentum distribution of 
the incident jet is broadened by collisions with scatter- 
ers. From the broadening of the transverse momentum 
distribution associated with the transverse momentum 
transfer due to the interaction of the jet with the ith 
scatterer, one infers that the corresponding longitudinal 
momentum transfer to the ith scatterer is given again by 
Qiz ^ l9iTp/2po as in Eq. © in the coherent collision 
limit. 

From the above review, we recognize that in the po- 
tential model, all recoils of the scatterers have been as- 
sumed to be zero initially to allow a static description of 
the centers of the scatterers. The scatterer longitudinal 
recoils are only subsequently corrected as an appendage 
to the dynamics of the transverse deflection of the inci- 
dent jet. Such a potential model may be adequate for 
incoherent collisions and for quasi-elastic coherent colli- 
sions in which the longitudinal momentum transfers to 
the scatterers are very small. 

As the longitudinal recoils of the scatterers are at the 
focus of our attention, it is important to realize that the 
potential model cannot be used to examine the longitudi- 
nal recoils of the scatterers in a general treatment where 
one wants to explore the behavior of the probability dis- 
tribution for the scatterer longitudinal momentum trans- 
fer over a large domain, or in regions where there may be 
a high probability for recoiling scatterers to share sub- 
stantial fractions of the longitudinal momentum of the 
incident jet. A general treatment of coherent collisions 
necessitates the use of the longitudinal recoils of the scat- 
terers as independent dynamical variables. The longitu- 
dinal recoils of the scatterers are however not allowed 
to be independent dynamical variables in the potential 
model. This leads us to forgo the potential model and 
to turn to the use of Feynman amplitudes for coherent 
collisions. 

In the Feynman amplitude approach such as given in 
what follows for a set of the initial momenta of the jet 
and scatterers, the (n -I- 1) final momenta involved in the 
collision process are dynamical variables. The Feynman 
amplitudes give the probability amplitudes for various 
reaction channels as a function of these dynamical vari- 
ables. In particular, the longitudinal momenta of the 
recoiling scatterers and the corresponding longitudinal 
momentum transfers from the jet to the scatterers are 
dynamical variables, in contrast to the potential model 
in which they are frozen and dependent quantities. Feyn- 
man amplitudes with explicit scatterer momenta on the 
external legs are the proper tools to examine the proba- 
bility amplitudes as a function of the longitudinal recoils 
of the scatterers. 

In the Feynman amplitude approach, the momenta of 
recoiling scatterers are intimately linked together with 
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the momentum of the incident jet in various Feynman 
ampUtudes and their corresponding symmetrized n\ per- 
mutations. The n scatterers and the jet comprise an 
{n + l)-body system, with 3(n + 1) degrees of freedom 
constrained by the four-dimensional energy and momen- 
tum conservation. Among the 3(n-|-l) degrees of freedom 
are the transverse and longitudinal momentum transfers 
from the jet to the scatterers, (qiTiQiz), * = 1, As 
the recoiling scatterers and the recoiling jet are tied to- 
gether by the coherent collision, they can share the ini- 
tial longitudinal momentum of the jet. We shall find 
from the Feynman amplitude approach that (qiz) can be 
a substantial fraction of the longitudinal momentum of 
the incident jet (see Section IV). The potential model 
in contrast gives the result qiz = \qfrp\/2po that cannot 
be relied upon for a general coherent collision, because 
(i) the potential model confines the system to regions of 
small Qiz , (ii) as it has been designed for elastic and quasi- 
elastic processes, the potential model does not allow q^z 
to be independent dynamical variables, and (iii) as a con- 
sequence, the potential model precludes the exploration 
into other qiz regions where the probability distribution 
as a function of qiz may be large in coherent collisions. 

In the general treatment of coherent collisions with the 
jet, we are therefore justified to forgo the potential model 
and use the Feynman amplitude approach to study the 
transverse and longitudinal recoils of parton scatterers. 



III. BOSE-EINSTEIN INTERFERENCE OF THE 
FEYNMAN AMPLITUDES 



(a) 



(b) 



p-q 



1 V 



l9i \92 



V P-ii |l 



1\ /^i 



p'=p-qfq2 

FIG. 1: Feynman diagrams for the collision of a fast fermion 
p with medium fermions ai, 02, with the emission and ab- 
sorption of virtual bosons of momenta q\ and 52. 



1(a) suffices. It can be cut into two disjoint pieces. The 
collision process consists effectively of a sequence of two 
two-body collisions. We shall show in Section IV that for 
a ladder diagram of the type shown in Fig. 1 (a) for quasi- 
elastic collision, with intermediate fast particles assumed 
to be on the mass shell and ^ 0, the (l+ri) body cross 
section is a product of n two-body cross sections. 

On the other hand, if the collision is characterized by 
A <C Azcoh the collision process is in a coherent (1+2)- 
body collision. There is an additional Feynman ampli- 
tude Mb for diagram 1(b), obtained by making a sym- 
metrized permutation of the bosons in diagram 1(a), 



-/u(p')7m; 
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We consider a jet p passing through a dense medium 
and making (coherent) multiple collisions with medium 
partons. The assembly of medium partons has an ini- 
tial momentum distribution. We choose to work in 
the center-of-momentum frame of the parton scatterers 
{oi, a2, a„}. We select the longitudinal z-axis to be 
along the momentum of the incident jet. 

As an illustration of the salient features of the interfer- 
ence of the Feynman amplitudes, we consider first mul- 
tiple collisions of a fermion p with rest mass m and two 
fermion scatterers oi and ai with rest masses mi and m2 
in p + ai + 02 p' + a'j^ + We study this problem in 
the Abelian gauge theory in this section and in the non- 
Abelian gauge theory in Section VI. We shall consider 
the collision of gluons in Section VII. 

Using the Feynman rules as in Ref. 18411 . the Feynman 
amplitude for diagram 1 (a) is given bv |95l. l96j 



Ma = -.9^1(^)7^ 



1 



77m"(p) 



X ^u{a2)'-f^u{a2)\u{ai')jf,u{ai). (9) 

qi qt 

If the spatial separation between the scatterers is so 
large such that A ^ A^coh the fast particle is nearly on 
the mass shell after the first collision and the diagram 



X — M(a2)7i/w(a2)— u(ai')7pu(ai). (10) 

The trajectories for diagram 1(a) and 1(b) arc both pos- 
sible paths in a coherent collision, leading from a set of 
initial states to a set of final states. By Bose-Einstein 
symmetry, the total amplitude M for coherent collisions 
is the symmetrized sum of Ma and Mi,. 

We consider the high-energy limit and assume the con- 
servation of helicity with 

{Po, \p\,Po, \P'\} > {\ar\,qo, kill > m, for i 1,2. 

We shall later find that {qiz)/Pz is of order l/2n. The 
above high-energy limit of {po, \p\,p'o, \p'\} > {go, 
provides a reasonable approximation to the gross features 
of the collision processes. Higher order corrections to this 
approximation will lead to refinements that will be the 
subjects of future investigations. 



In this high-energy limit, we have approximately 84 1 



u{a )7m^(q)-W ,^ + \ \ = — (H) 



-luu{p) 



2p ■ qi — ie 

"^PuPn 



■7.w(p),(12) 
(13) 
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where e is a small positive quantity. We shall be inter- 
ested in the case in which the ferniion p' after the collision 
is outside the medium and is on the mass shell. The mass 
shell condition can be expressed as 



{p- qi- q-iY - "^^ 



-2p-qi-2p-q2 - 0. (14) 



The symmetrized sum of the Feynman amplitudes Ma 
and Mf, in the high-energy limit is 



Mr 



1 



1 



g'^ 2p ■ di 2p ■ 0.2 

2m miq\ m2(7| \2p ■ qi — ie 2p ■ q2 — ie 



(15) 



Note that the amplitudes Ma and Mf, correlate with 
each other because of the mass-shell condition ([T4)) . The 
real parts of the amplitudes destructively cancel, and the 
imaginary parts interfere and add constructively, to re- 
sult in sharp distributions at p ■ qi ~ and p ■ q2 ^ 0, 



M 



g'^ 2p ■ di 2p ■ 0,2 ( . 



2m miqf 77i29| 
where 



inA{2p-qi) + iTTA{2p-q2)y (16) 



A{2p-qi) = -- 



7r(2p-(7i)2 + e2' 

which approaches the Dirac delta function 6{2p ■ qi) in 
the limit e ^ 0. 

We consider next the case of the (coherent) multiple 
collisions of a fast fermion with three fermion scatterers 
in the reaction p + Oa + 02 + p' + a[+ a'2 + Og. The 
Feynman diagrams for the multiple collisions process are 
shown in Fig. 2. 
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{2p- 


92 


- ie){2p ■ q2 4 
1 


~2p- 


93 - 


-^e) 


{2p 


• 93 


- »e)(2p • 92 - 
1 


f 2p 


■ 93 


-^e) 


{2p- 


93 


- »e)(2p ■ 93 -\ 
1 


-2p- 


91 - 


-^e) 


{2p 


■ 91 


-ie)(2p-93- 


f 2p 


■ 91 





(17) 



Each amplitude term in the curly bracket contains a dis- 
tribution that is non-zero in various regions of p-qi . These 
Feynman amplitudes interfere among themselves. We 
can write Eq. ()17|1 in the form 



M 




2m I ^ ^ m.o; 
U=i ' 



j=i ELi 2p • 9^ - 
symmetric permutations >. (18) 



Generalizing to the case of the coherent collision of a 
fast fermion with n fermion scatterers in the process 

p + oi + ... + a„ — > p' + a[ + ... + a^, (19) 

the total Feynman amplitude is 

i'^2p-a,'i 



M 



9 



n 



2m TJiiqf 



■■A^(9i,92,---,9«), 



(20) 



where A^(qi, 52, 9n) is the sum of n! amplitudes in- 
volving symmetric permutations of the exchanged bosons 
given by 



7W(gi,g2,---,9n) 



n 



1 



t=i ELi 2p • 9» - »e 
symmetric permutations. (21) 



The above sum involves extensive cancellations. Remark- 
ably, it can be shown that this sum of n\ symmetric per- 
mutations turns out to be a product of sharp distribu- 
tions centered at 2p • ~ [8J| , 



FIG. 2: Feynman diagrams for the multiple collisions of a fast 
fermion with three medium fermions in different permutations 
of the exchanged bosons. 

In the high-energy limit, the Feynman amplitude for 
the collision is 



M 



g^ 2p ■ di2p ■ 02 2p ■ 03 1 
2m mi 7712 ma 9i9293 
1 

{2p ■ qi ~ ie){2p ■ qi + 2p ■ q2 - ie) 



X(gi, 92, .., 9«)A(^2p • q,) ^ (2^i)"-i[]A(2p • q,), (22) 

i=l i=l 

which gives 

Al((7i,g2,...,9n) = ^^^ Y\ n^(2p.g,)j. (23) 

1=1 ^j = lj-5^i 

Eqs. ([2T|) - ((23| indicate that for coherent collisions, there 
is a destructive interference of the Feynman amplitudes, 
leading to a cancellation in most regions of the phase 
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space. However there is a constructive interference in 
some other regions of the the phase space, leading to 
sharp distributions at 2p ■ qi 0. The constraint 
A(2p ■ Qi) is mathematicany the same as the constrain 
A((p — Qi) — m?) in the high-energy Hmit, which can be 
interpreted as the intermediate state pi+i={p — qi) to be 
on the mass sheU. It should however be kept in mind that 
the intermediate state Pi+i = {p — qi) in the constraint 
of Eqs. (gni), (1221), and is not a single particle but 
a quasi-particlcs that is the result of the interference of 
many amplitudes. 

Equations ((20)) - (|23)) for the total Feynman amplitude 
for a coherent collision as a product of delta functions of 
2p-qi are similar to previous results obtained for the emis- 
sion of many real photons or gluons in bremsstrahlung, 
and for the sum of ladder and cross-ladder amplitudes in 
the collision of two fcrmions I84i-l90ll. 



IV. CROSS SECTION FOR (l+n)-BODY 
COLLISION 



We wish to obtain the cross section for the coherent col- 
lision of a fast fermion p on n fermions {ai, 02, 03, a„} 
in reaction While the Feynman rules to construct 

the Feynman amplitude are well known, the differential 
cross section for the reaction ((T9)) as a single (l-l-n)-body 
collision process, written in terms of the Feynman ampli- 
tude M, does not seem to be available in the literature. 
We would like to examine the cross section rules in some 
details in this section. 

It should be mentioned that the cross section rules for 
the (l-l-n)-body collision depend on the Feynman rules 
for the Feynman amplitude M. We shall use the Feyn- 
man rules of Refs. [SJ, for M, which adopt the nor- 
malization uu — 1 for the fermion wave function. A 
different normalization of the fermion wave function ac- 
cording to uu = 2m, as given in Ref. [93|, will lead to 
slightly different cross section rules. 

To be concrete, we use the auxiliary variable pi with 
i = l,...,n to label the intermediate state (or quasi- 
particle) of the incident fermion prior to its exchange 
of a virtual boson with scatterer ai (or after its exchange 
of a virtual gluon with scatterer ai_i), with the special 
case pi = p and Pn+i = p' ■ The differential cross section 
for the collision process ([T5| involving (l+n) fcrmions is 
given by 



the scatterer of rest mass 



dV2m Dp'jp') f -A- d^a[2mi Dija',) 
(27r)3 2p'o \l\ (2^)3 2a,;o 



(24) 



where M is the total Feynman amp litude as follows from 
the Feynman rules of Refs. [HjIqI^H^. The quantity fpi 
is the flux factor for a two-body collision between pi and 



fpi 



4^ {pi ■ aiY - (to m~Tf 
(2to)(2to^) 



(25) 



and rriiT = \/ mj 



''iT- 



In Eq. (|24p we construct the flux factor for the (1-l-n)- 
body collision by requiring that it should be the same for 
all scatterers. Hence, it is the product of flux factors from 
all scatterers, n"=i fpi- We also need to require that in 
this flux factor, the momentum of the incident fermion 
and the scatterer fermions appear in the same power. 
Thus, the flux factor need to be linear in the energy of 
the colliding particles. But the product of nr=i /p* gives 
nr=i?'*o- It is therefore necessary to multiply 0"=! fpi 
by the factor Y['i=2i^/Pio) get the flux factor for the 
(l-l-n)-body collision that is linear in the energy of all 
colliding particles, resulting in the factors as given in the 
cross section formula in Eq. ([211). 

The quantity T^, with {i = 2, ..n} is the collision time 
or the mean lifetime of the intermediate state (or quasi- 
particle) pi of the incident fast particle prior to the inci- 
dent particle exchanging a virtual boson with (or after 
the incident particle has already exchanged of a virtual 
gluon with scatterer a^-i). It arises in the cross section 
formula because the cross section is proportional to the 
transition rate per unit time of collision and thus de- 
pends on the lifetime of the intermediate state pi of the 
projectile before it exchanges a virtual boson with the 
scatterer a^. The precise value of Ti needs not concern 
us at this point as it will be canceled later by the delta 
function that specifles the energy-momentum constraint 
of the intermediate state pi, evaluated at the condition 
of the constraint of the delta function. 

The states of the incident particle p' and the medium 
scatterer i' after the collision can come in different forms. 
They can be described by the functions Dpr /2pQ and 
Di/2a^Q. For example, the final fast particle p' resides 
outside the medium after the collisions without interac- 
tions. Its state can be adequately described as being 
nearly on the mass shell. 



Dp'ip') S{p'o - Vp'^ + m2) 



'^P'o 



S{p'^ - to2). (26) 



It may subsequently break up as jet fragments in a col- 
limated narrow cone, in accordance with a jet fragmen- 
tation function. We shall examine the states of the scat- 
terers in the next section. 

To show the validity of Eq. (p4l) . we note first that 
it is clearly correct for the case of a two-body collision 
with n = 1, as given in Appendix A-3 of Ref. [9^. We 
would like to test whether it gives the correct result for 
the case of quasi-elastic collisions where the answer is 
known. Such of quasi-elastic collisions is described 

by the generalization of ladder diagrams of Fig. 1(a), and 
2(a), with the incident particle nearly on the mass shell 
and qio 0. The diagram can be effectively cut at the 
intermediate fermion lines, resulting in a sequence of n 
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disjoint two-body collision diagrams, and the many-body 
cross section is just the product of n two-body collisions. 

We would like to show that for the above quasi-elastic 
collision case the Feynman amplitude approach with the 
cross section rules of Eq. (|24p also leads to the known 
result. 

In the high energy limit, the Feynman amplitude for 
the ladder diagrams of Fig. 1(a), and 2(a) can be gener- 
alized for n scatterers as given by 



where the mean lifetime T^+i cancels the corresponding 
factor in the denominator of Eq. (|30p. Using Eqs. (PT|) 
and ([32]) in Eq. dSO]), we obtain 



^ n^^(2rn)27rA(p?+,-m2)-'^''^^2"^^^^^«^) 



1=1 



pi 



fpn 



i=l 



2m 



'n-l 

n 



2pi ■ hi 



1 



mnql 



(27) 



rfV2m Dp,ip') d4a;2m„ D„{a'J 
(2^)3 2p'o (27r)3 2a„o 



(27r)3 2a,o 



(33) 



For quasi-clastic colhsions with q,o - 0, the incident par- Inserting the unity factor JJ"^^ 5{p, + a, p,+i - 
tide is nearly on the mass shell with ^ m^. As "■'i)'^ -1 into the above equation and noting that 



a consequence, the propagator represented by l/{pf_^i — 
■m? + ie) in the Feynman amplitude can be approximated 
by the pole term 



The Feynman amplitude becomes 



M 



g j 

2^ I n 



2pi • di . 



i7rA(pj — m ) 



2p' ■ CLn 

mnql 



(28) 



(29) 



We substitute this Feynman amplitude in Eq. ([24| . and 
we obtain 



(To 



1 



9 



An 



{OLi /pj{nr=2 (™/Pio)rj (2m; 



Tl-l 



2 

'4„' 



1=1 



? 4 



d-^a^^2m, A«) 
(27r)3 



2a 



iO 



(2p ■ Qn 

2 4 



-{27:r5Hp + J2a,-p'-Y,a^) 



dV2m jJp>(y) d4a;2m„ i^„(a;,) 
(2^)3 2p^ (2^)3 2a„o 

On the other hand, for the two-body collision pi + ai 
Pi+i + a'i, the Feynman amplitude is 



(30) 



the differential cross section for pi + ai ^ Pi+i + is 

^2B|2 



da: 



2B 



-{2tt) 6 [pi + fli - pi+i - a'i) 



fp 

d'^Pt+i2m A+i(Pt+i) d'^a'^2mi Di{a'i) , , 



(2^)3 2p(,+i)o (2^)3 2a. 



iO 



we can rewrite Eq. p3p as 



1=1 



2B 

pi 1 



(35) 



which is the cross section for quasi-elastic collisions. This 
show that Eq. (|24p gives the correct cross section for the 
case for which the answer is known. It can therefore be 
used for a general treatment for the cross section in terms 
of many-body Feynman amplitudes in the general case of 
coherent collisions. 

Equation ([M)) exhibits explicitly the rules for the cross 
section for coherent collision of an incident fast fermion 
with n fermion scatterers. The other cases for incident 
gluons and gluon scatterers will be taken up in Section 
IX. 



V. CONSEQUENCES OF THE BE 
INTERFERENCE ON THE RECOIL OF 
FERMION SCATTERERS 



M. 



2B 



2pi ■ a» 
2m rriiq^ 



(31) 



Upon writing A(p|_|_j^ — m?) in terms of the mean lifetime 
Ti+i = l/Fi+i of the intermediate state oi pi+i, we have 



A(p-+i - rr?) 



we get 



ri/(47rp(j+i)o) 



P(i+1)0 



Pi+i + 



+ rf+i/4 



|A(p?^i-m,2)P = A(p?^i-m2) 



2T,: 



i+l 



, (32) 



The states of a medium scatterer after the collision can 
come in different forms. As it acquires a large amount of 
energy and momentum from the incident jet, the medium 
scatterer can be far off the mass shell and it may subse- 
quently fragment into a cluster of on-mass-shell particles. 
For an incident jet of energy 10 GeV with scatterers in 
a medium, it is however unlikely that the final scatterers 
themselves be energetic enough or far off the mass shell 
to fragment as clusters. Far more likely is the case of 
the scatterer a^ to acquire a fraction of the incident jet 
energy and momentum and become only slightly off the 
mass shell, with the degree of its being off-mass shell de- 
scribed by a width F^ . As the medium particles reside in 
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an interacting medium, it will be subject to the interac- 
tions of the medium before and after the collision. We 
can describe the mass-shell constraint for the of medium 
scatterer a' after the collision as 



[a-o - ^/{a',-g^A)^ + im, + S)^ + g,Ao]^ + -f 



where A = {A, Aq\ and S are the vector and scalar mean 
fields experienced by the medium scatterer after the 
collision, respectively. As the mean fields and scatterer 
widths increase with medium density and are presumably 
quite large and dominant for a dense medium, we shall 
approximately represent Di as an average constant that is 
only a weak function of a[. Other descriptions of Di/2a[^^ 
for the states of the scatterers are also possible but may 
not be as general; they can be the subjects for future 
investigations. 

We integrate over 5{p' + ^ qi —p)d^p' and change vari- 
ables from a[ to qt = a[ — Ui. Writing out the matrix M 
explicitly as in Eq. ^U\\ and using the relations in Eq. 
and dial), we obtain from Eqs. pi]) 



'(2m)'^ 



(27r)"+i 



:/:.{2p-q^)dq^,\\[—^ 



.4=1 



(36) 



where for simplicity, we have taken = m. The distri- 
bution A(p • qi) can be written as 



A(2p • q, 
where C is 



1 



PO+Pz 



C 



-(Po - Pz){qm + qiz) + 2pr ■ 

PO+Pz 



(37) 



(38) 



and the quantity ^ approaches zero in the high-energy 
limit of large po- The function A(gio ~ Qiz — C) provides 
the constraint at qio — qi^ ^0. As a consequence of the 
constraint, the integration A{qiQ — qi^ — QdqiQ can be 
carried out, yielding 



A(2p • qt)dqtO ' 



^ A{qio-qiz + OdqiQ = ^. (39) 

Pz +Po ^Pz 

Noting that A(2p • qi)^{p-qif-m'^ , we can express 
A(2p- qi) in terms of the mean lifetime T^+i = l/Ti^i of 
the intermediate state p — qi= Pi+i of the fast projectile, 
after it has exchanged of a virtual boson with the scatter 
I, as 



A(2p • q{) = 



ri+i/47rp( 



i+l)0 



[P( 



i+l)0 



then we have 



|A(2p.g,)l'- A(2p.9,)- 



T^Plyi+V)^ 



(40) 



(41) 



The above mean lifetime T^j^x cancels the mean lifetime 
Ti in the denominator of Eq. (j36p . The boson propagator 
g| in the denominator becomes 



(li = (QiO + qiz){.qio - Qiz) - kirP ~ -|9»tI 
We obtain 



If \ " f -A- 8A (2p • aj)2 dqi^dq^rj.] 



4 \mp2 



n 



(42) 



(43) 



where a=g^ /A-k. The last factor in the curly bracket is 
dimensionless, and d"(T has the dimension of {o? jmpzY^ ^ 
as it should be. 

To find the probability distribution for the longitudi- 
nal momentum transfer qi^., we introduce the fractional 
longitudinal momentum kick 



Xi = —, dq.,^ = Pzdxi 

Pz 



(44) 



To investigate the Xi dependence of the factor 
(2p-ai)2/2aJo in Eq. (gSl), we note from Eq. dTT]) that 
hi can be written as a function of q^ and a^. 



a- ^ / "io + "'^ g' I g^o + m + ato + TO aj ^^^^ 
V a'm + m 2 ^(a'^g + m)(a,o + m) 2 ' 

Because of the A(2p • qi) constraint, the factor 
(2p-ai)2/2a^g in Eq. (gSl) becomes 

(2p ■ a,)2 ^ (a;.Q -f flio)^ (p • a,)2 _ (p-a»)^ ^^g^ 



2a: 



iO 



2(a:;o)2 



We obtain from Eq. (gSI 



8 A Ki(p • Oi)^ 




, . , J I] ITlQiO 

dxidx2---dxndqirpdq2T---dqnT 



\QnT\ 



(47) 



The fermion scatterers can possess different initial ener- 
gies a^o at the moment of their collisions with the en- 
ergetic jet. In the case when a^g ^ (?iOi the factor 
approaches 1/2 -f 0{aiQ/qiz) with (aio/qiz) ^ 1. In the 
other extreme when a^o ^ ftoj the factor k approaches 
2 + 0{qiz/aio) with {qiz/ciio) <C 1. The dependence of Ki 
on and qjjn is weak in either limits, and such depen- 
dencies can be neglected in our approximate estimate. 
We obtain then 



d"(T - C, 



where Cq is 



dxidx2...dxndqirpdq2T---dqnT 



Ca 




(48) 



(49) 
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the quantity Cq is a weak function of Xi and q^j, and can 
be approximated as a constant. 

By symmetry, the fraction of momentum transfers Xi 
for different scatterers should be approximately the same 
on the average, and x™^'^^ 1/n. Then as far as Xi is 
concerned, the average distribution is 



dP 1 

n e( x^), 

axi n 

and the average longitudinal momentum fraction is 
{X.) - ^ or (g,.) ~ -. 



(50) 



(51) 



The above result in Eq. ([50| indicate that the probabil- 
ity distribution of the longitudinal momentum transfer of 
a scatterer is approximately flat in Xj. The quasi-elastic 
scattering with qiz occurs, but with about the same 
probability as other longitudinal momentum transfer up 
to qiz = Pz/n. On the average, a scatterer acquires ap- 
proximately l/2n fraction of the incident jet longitudinal 
momentum. 

Because of the recoils of the scatterers, the incident jet 
loses its longitudinal momentum. On the average, the 
jet loses a longitudinal momentum fraction of l/2n after 
each collision with a scatterer. When the jet emerges 
out of the medium after colliding with n scatterers, the 
(average) fractional jet longitudinal momentum loss due 
to scatterers recoils is approximately 



that the fermion scatterers emerge in the direction of the 
incident particle, each carrying a fraction of the forward 
longitudinal momentum of the incident particle that is 
inversely proportional to twice the number of scatterers, 
{qiz) ^ Pz/'2^n. 

We have presented an explicit derivation of the differ- 
ential cross section as a function of the momenta of the re- 
coil scatterers by making many simplifying assumptions. 
The explicit derivation has the advantage that it allows 
future refinements on some parts of the calculation by 
modifying some of the simplifying assumptions. 



VI. BOSE-EINSTEIN INTERFERENCE FOR 
COLLISIONS IN NON-ABELIAN THEORY 

The above considerations for the Abelian theory can 
be extended to the non- Abelian theory. As an example, 
we consider a quark jet p making coherent collisions with 
quarks oi and 02 in the reaction p + oi + 02 — > + a'j^ + 
a'j, in the non- Abelian theory. We shall neglect four- 
particle vertices and loops, which are of higher-orders. 
The Feynman diagrams are then the same as those in 
Fig. 1. One associates each quark vertex with a color 
matrix T^^^'"^^ where the superscript p, 1, or 2 identifies 
the quark p, oi, or 02, and the subscripts a or /3 give the 
SU{i) color matrix index. The Feynman amplitude Ma 
for diagram 1(a) is 



(fractional jet longitudinal momentum loss) 



2n 
(52) 



Having obtained the differential cross section Eq. 
we can infer the distribution of the scatterers with respect 
to the incident particle axis. Equation pS)) indicates that 
the reaction has a high probability for the occurrence 
of small values of Iqi^l, in the passage of an energetic 
fermion making coherent collisions with medium partons. 
The singularities at jq^yl ^ in Eq. ([iS]) correspond to 
the case of infrared instabilities that may be renormal- 
ized, and a momentum cut-off Acut may be introduced. 
As q.^rp are independent degrees of freedom, Eq. (|43|) or 
PS)) shows that the standard deviation of the transverse 
momentum distribution of the scattered incident particle 
is related to the transverse momentum transfers to the 
scatterers by 



((Pt)') = ((E9^t)')=E' 



n{qiT) 



(53) 



as in a random walk in the transverse direction. 

Equations P5)) - ([5T|) indicate further that the scatterers 
acquire an average longitudinal momentum {qiz) ^ Pz/'in 
that is expected to be much greater than (|qir|). Thus, in 
a coherent collision there is a collective quantum many- 
body effect arising from Bose-Einstein interference such 



1 



I — — m + 26 



X ^u(a2)T^^^7i.M(a2)^u(ai')r^^^7M"(<^i)- 
The Feynman amplitude Mb for diagram 1(b) is 



1 



— m + le 



2 



In the high-energy limit for coherent collisions, the sum 
of the Feynman amplitudes is 



M 



g 2p • oi 2p • 02 1 
m q^qf 



2m m 



' rp(p)rp{p)rj.{2)rp(l) 

2p ■ qi — ie 



rp(p)rp{p)rp{2)rp{l)^ 

J-a J-^ J-^ J-a 

2p- q2 - ie 



.(54) 



We can rewrite the product of the color matrices for the 
quark jet p as 



(55) 



tL^^T^/^ = I {[Tlr\TiP^U - [Tlf\m-) ■ (56) 



11 



We label the propagators in Eq. as A^^ and Aib, 

(57) 



Ma = -, and Mb = ^ 



2p ■ qi — ie 2p ■ q2 ~ ie 

The Feynman amplitude for coherent collision is then 
2p ■ ai2p ■ 0.2 1 



M 



2m m m q2li 
xliMa+Mb) 



+ {Ma-Mb) 



[l p ,1a l-lp la 



(58) 



The first term inside the curly bracket has the same 
space-time structure as what one obtains in the Abelian 
theory. It is given by the Abelian Feynman amplitude in 
Section 11, multiplied by the color factor 



rpi2)rp{l) 

2 ^ P 



(59) 



The sum of A^^ and Mb leads to sharp distributions at 
p • (7i ~ and p • 92 ~ 0, 



Ma + Mb^ iA{2p ■ qi) + iA{2p ■ ga)- 



(60) 



The second term is new and occurs only in the non- 

Abelian theory, as it involves the commutator of Tj^^'' 
iv) 

and Ta ■ It involves the difference oi Ma and Mb in 
which the sharp distributions cancel each other, leaving 
a broad distribution, 



Ma -Mb 



Ap ■ qi 



(2p.gi)2 + e2- 



(61) 



From the above analysis, wc find that the color degrees 
of freedom in QCD bring in additional properties to the 
Feynman amplitude. Bosc-Einstcin symmetry with re- 
spect to the interchange of gluons in QCD involves not 
only the space-time exchange symmetry but also color 
index exchange symmetry. The total exchange symme- 
try can be attained with symmetric space-time ampli- 
tudes and symmetric color index factors as in the first 
Ma + Mb term in Eq. ([55|. The total symmetry can 
also be attained with space-time antisymmetry and color 
index antisymmetry, as in the second Ma — Mb term in 
Eq. dSSD . 

We consider next the case for the collision of a quark 
jet with three quark scatterers in the reaction p + aa + 
02 + 03 — > p' + a[ +03 + 03 in the non- Abelian theory. 
The Feynman diagrams for the collision process are the 
same as those in Fig. 2 where we associate the color ma- 
trices Ta^\^ ^i^'; Tc^^ of color indices a, b, c with fcrmion 
scatterers oi, 02 and 03 respectively. In the high-energy 
limit, the Feynman amplitude for the collision is 



AI = 



2p ■ ai2p ■ a22p ■ 03 
2m "^"^ 91 92 93 



TaTbTcMi23 + TaTcTbMi32 + TbT,TaM231 



+TbTaT,M213 + T,TaTbM3i2 + T,TbTaM321^^2) 

where Ta^b,c without a superscript are the color matri- 
ces for the incident quark jet p, and the amplitudes 
Mi23,Mi32,M23i,M2i3,M3i2,M32i are sequentially the 
six terms in the curly brackets of the Abelian Feynman 
amplitudes in Eq. (|17p . It is easy to show that the quan- 
tity in the curly bracket of Eq. (p^ can be re-written 
as 

M = {Ta[Tb,T,]++Tb[T,,Ta]++T,[Ta,Tb] + } 

X [M123 + Mi32 + Af231 + M2I3 + M3I2 + A/32l]/2 

- Ta[Tb, Te]+[M231 + Af213 + ^312 + M32l]/2 

- Tfc[r„ Tj + [Mi23 + Afi32 + A/312 + A/32l]/2 

- Tc[Ta,Tb] + [Mi23 + A/132 + A/231 + A/213] /2 
+ T,[rb,Tj_ (A/123- A/132) /2 

+ Tfc[r„ TJ_ (A/231 - A/213) /2 

+ T,[r„rb]_(M3i2- A/321) /2. (63) 

The first term on the right hand side involves the sym- 
metric sum of the space-time part of permuted ampli- 
tudes as in the Abelian case, multiplied by the symmet- 
ric permutation of the color indices. It yields a Feynman 
amplitude that is just the Abelian Feynman amplitude 
multiplied by the color factor 



7.(1)7^(2)7^(3) 

r< lb -'c J 71 TT^ 7. " 
'^r,^= Z <la[lb,lcl+- 



-Tb[T„Ta]++T,[Ta,Tb]- 



The other terms involve partial symmetry and antisym- 
metry with respect to the exchange of color indices. Sim- 
ilar studies on the collision of a jet with n partons can 
be carried out as in (ssl - fooj . For our present work, it 
suffices to note that there will always be a component 
of the Feynman amplitude that is symmetric under both 
space-time exchange and color index exchange involving 
the sum of all space-time amplitude components, similar 
to the Ma+ Mb sum in Eq. ([55]) and the first term on 
the right-hand side of Eq. . There will also be other 
space-time antisymmetric and color index exchange an- 
tisymmetric components. 

For the space-time symmetric and color index exchange 
symmetric component, the Feynman amplitude is equal 
to the Abelian Feynman amplitude multiplied by a color 
factor. It will exhibit the same degree of Bose-Einstein 
interference as in the Abelian theory. Previous analy- 
sis on the longitudinal momentum transfer of recoiling 
fcrmions in the Abelian theory in Section 2 can be ap- 
plied for the non- Abelian theory for this space-time sym- 
metric and color index exchange symmetric component. 
There is thus a finite probability for the presence of delta 
function constraints to lead to recoiling quarks receiving 
significant moment kicks along the direction of the inci- 
dent quark jet. 
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VII. COLLISION OF A GLUON JET WITH 
QUARK SCATTERERS 

It is of interest to generalize the above considerations 
to the coherent multiple collisions of a gluon jet. We 
shall neglect four-particle vertices and loops, which are 
of higher-orders. The Feynman diagrams for the collision 
of a gluon jet with medium quarks or medium gluons then 
have structures and momentum flows the same as those 
in the collision of a quark jet with quark scatterers. In the 
high-energy limit, the propagators and the three-particle 
vertices have approximately the same momentum depen- 
dencies, and the Bosc-Einstein symmetry with respect to 
the interchange of the virtual bosons is the same. One 
expects that aside from the presence of color factors and 
color indices, the results for the Bose-Einstein interfer- 
ence in collisions with a gluon jet or a quark jet should 
be similar. This is so because the high-energy processes 
are insensitive to the spins of the colliding particles, as 
the current carried by a high energy particle is dominated 
by its center-of-mass motion, much more so than its spin 
current [H, [oi, [gi] . 

It is instructive to study the coherent collision of a 
fast gluon p with two quarks in non-Abelian gauge field 
theories in the reaction p + ai + 02 — > p' + a'j^ + a'2 as an 
examples of the type of BE interference for the collision 
of a gluon jet. The collision process is represented by the 
two Feynman diagrams in Fig. 3. 



p 


p-qi 


P' 


aa • 




> C7 




In \ 








>ee 


ai ] 


> ffl'l 12 I 


> "2 




hrj kn 

(a) 


IX 




xu{a'2)^^u{a2)u{ai')^su{ai), (64) 

where e{p) and e(p') are the polarization vectors 
for gluons p or p' respectively. In the high-energy 
limit in which |p| |aj|, m}, terms of order 

{|'Zil/|pM'3'i|/|p|j ill the three-gluon vertices can 
be neglected, and the helicities can be assumed to con- 
serve. We then obtain 

{p' , a'lhrj ^ a^lX \ Ma \ paa,aif(j),a2kK) 



TO 



(65) 



where the coefficients fbec and fbad in Eq. (j64p have been 
expressed as matrix elements of matrices Tjf'' and tJ''-' 
of the incident gluon p (or p') between gluon color states 

{T^%a = ifdba, (66) 

(TiP))c6 = ifecb. (67) 

Equation (|65p can be re-written in a matrix form as 



Ma^g 



,2p-a^2p-~a2T^^^Tl^^TPTl^^ 



9 9 2 



2p- qi 



(68) 



We can obtain a similar result for the Feynman ampli- 
tude for diagram 3(b) by permuting the vertices of the 
exchange bosons. As a consequence, the sum of the two 
Feynman amplitudes from Figs. 3(a) and 3(b) is given by 



fcj) hrj k.K IX 
ib) 



FIG. 3: Feynman diagrams for the collision of a gluon jet 
p with quark partons oi, 12 with the momentum transfers 
qi and q2 in the reaction p + ai + ai — >■ p' + a'l + 02. Here, 
{a, b, c, k, 1} label the color states and {a, /3, 77, A} label 
the momentum components of various particles. 

We consider the Feynman amplitude matrix element 
between the initial three particle state \paa, aif4i^ a2kK) 
and the final state \p'cj, a[hri, a^lX) where {a, /, fc, c, /i, V\ 
label the color states and {a, (j), k, 7, 77, A} label compo- 
nents of various particle currents. Using the Feynman 
rules and the overall phase factors as given in Ref. [s^ . 
the Feynman amplitude matrix element of Fig. 3(a) is 
given by 

{p'cj,a'ihri,a'2lX \ Ma \ paa,aif(f>,a2kK) 

4 fbecfbad /rj^(2)\ / r\ 

x{ff7/3(p' + P - + g,Be{-p + qi - q2)-y + gcfiq2 -p')i3 

x{5/3a(2p - qi)5 + ga6{-p - qi)l3 + g5p{2qi - p)a}f*a{p) 



2p ■ di2p ■ 0,2 

in^qlql 

-'e 

2p ■ qi — ie 



rpip)rp(p)rp{2)rp{l) 

^d ^d 

2p ■ q2 - ie 



(69) 



We note that the above equation for the collision of a 
gluon jet is in the same form as Eq. (|54p for the collision 
of a quark jet except with the modification that in the 
above equation for a gluon jet, the operator T^p^ has 
matrix elements between gluon color states whereas the 
operator T^^*^ in Eq. ([54| for the quark jet has matrix 
elements between quark color states. As in Eq. ([M]) . we 
can likewise express the product of the color matrices as 



The Feynman amplitude for a gluon jet is then 
^2p ■ ai2p ■ 0.2 



(70) 
(71) 



M 



Oil 

m^qfqi 



{Ma+Mb)^^ 



(P) rp(P)^ rj.(2)rr.{l) 



} 



+ {Ma-Mb) 



(72) 
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Therefore, in the colhsion of both a gluon or a quark 
jet with quarks, there will always be a component of the 
Feynman amplitude that is symmetric under both space- 
time exchange and color index exchange, involving the 
sum of all space-time amplitude components. The simi- 
larity is so close that previous results concerning a quark 
jet in collision with quark scatterers apply equally well 
to a gluon jet. 




VIII. COLLISION OF A GLUON JET WITH 
GLUON SCATTERERS 

There is however a small difference in the coherent col- 
lisions of a jet with gluon scatterers. We can consider the 
collision of a fast gluon p with two gluon scatterers oi and 
02 as shown in Fig. 4. 

The matrix element of the Feynman amplitude in Fig. 
4(a) is given by 



FIG. 4: Feynman diagrams for the interaction of a gluon 
jet p with gluon scatterers ai, a2 with the momentum trans- 
fers gi and q2. Here, {a,b,c, ...,k,l} label the color states 
and {q, /3, 77, A} label the momentum components of vari- 
ous particles. 



{p' c^,a'ihrj,a'2l\ \ Ma 

3 f bee f bad f kief fhd 

qlq^C^p ■ qi - ie) 

x{ff7/3(p'+P-9l)e +9fle{- 
^{gpai'^P - qi)5 + 9a5{-p 

y-{gK.x{-a2 - a'2)c + gxM'2 



paa, aif4>, a2kK) 



the same form of the color operators, with the color oper- 
ators that operate on the quark color states switched to 
operate on the gluon color states, and (3) change the mo- 
mentum variable di of Eq. (jlip for a quark scatterer to di 
of Eq. (|75p for a gluon scatterer when a quark scatterer 
is replaced by a gluon scatterer. The above Rule (1) and 
e^(p')e* (p)eA(a2)e* (a2)e^(a'^)e^(ai) (2) are well known results that were obtained previously 

on page 198 of Ref. [sJl- Rule (3) is for our case in the 
92)7 + 567(92 — p')/3} collision with quark or gluon scatterers. 
gsp(2qi — p)a} ^'^^ completeness, we can use these rules to obtain the 

/ ^ ■) I Feynman amplitude for the collision of a quark jet with 
gth[ <Z2 0'2)x} gluon scatterers as 

gs,t>{-qi +ai)^}.(73) 



p + qi 

- qi)i3 ^ 
+ g2)K - 
+ 91)0 



In the high-energy limit in which \p\ ^ {I^J, ™} and 
the helicities can be assumed to conserve, we obtain 

{p'c'y,a[hr],a'2lX \ Ma \ paa,aif(l),a2kK) 

(T^(P)\ l^(P)\ /'T^(2)\ IrpWx 



M 



Ap ■ aiAp ■ Q2 
2m qlqj 

2p ■ qi — ie 



rpip)rp{p)rpi2)rp{l) 

^d J-e ±d 

2p • 92 - 



.(77) 



= g 



qlql{2p ■ qi - ie) 



Ap ■ aiAp ■ a2 



where 



(74) 



(75) 



We can obtain a similar result for the Feynman ampli- 
tude for diagram 4(b) by permuting the vertices of the 
exchange bosons. As a consequence, the sum of the two 
Feynman amplitudes from Figs. 4(a) and 4(b) is given by 



M 



^Ap ■ aiAp ■ a2 



IX. LONGITUDINAL MOMENTUM 
TRANSFER FOR GLUONS SCATTERERS 

The cross section rules for Eq. ([5(7| has been written 
for the collision of involving fermions (or quarks). The 
cross section rule involving incident gluons and/or gluon 
scatterers can be generalized by the rule that in replacing 
a fermion with a gluon, the 2TOformion factor for a fermion 
is replaced by unity for a gluon. Thus, we have explicitly 
the cross section formula for the coherent collision of an 
incident gluon with n fermion scatterers given by 



g 



2 9 

qiqi 



9+nf 



y(p)y(p)2^(2)y(i) 



2p-qi 



+ 



rpiP)rpip)rpi2)rpW 



2p-q2 



(76) 



{n:Li/pj{n.:i2(i/2Ko)}{nL2?^.} 



By comparing the results for a quark jet on quark scat- 
terers [Eq. (|54p ] with those for a gluon jet on quark scat- 
terers [Eq. I69|)] or on gluon scatterers [Eq. ([76])]. we ob- 
tain the following simple rules to obtain from the results 
of quark jet on quark scatterers of Eq. ([54| to those in- 
volving gluons: (1) provide an overall multiplicative fac- 
tor of 2m when a quark is replaced by a gluon, (2) keep 



(27r)3 2p'^ (27r) 



d'^a[2mi Di{a[) 



2a 



iO 



where 



/pi (gluon on fermion) = 



4p ■ aj 
2m,- 



(78) 



(79) 



Similarly, the cross section formula for the coherent colli- 
sion of an incident gluon with n gluon scatterers is given 



14 



by 



W U g-\.jig 



where 



{nr=i/pj{nr=2(i/2Ko)}{n:L2 7^j 



(81) 



(27r)3 2p'^ lfi(27r)3 2a,o J' 



/pi(gluon on gluon) 4p • a,;. 



The cross section formula for the coherent colhsion of an 
incident quark of rest mass m with n gluon scatterers is 
given by 



g+ng 



{nr=i /pj{n.:i2 ("VKo)}{nr=2 



rfV Dp-ip')2m.f-Pj d^a[ D,{a[) 



(27r)3 2p^, l,i-^(27r)3 2a,o 



n 



where 



/pi(fermion on gluon) = 



4p ■ Uj 
2m 



(82) 



(83) 



The results in the last section indicate that for the 
Feynman amplitude there is a modification from 5 of Eq. 
(|TT]) for a quark scatterer to of Eq. (|75p for a gluon 
scatterer, when a quark scatterer is replaced by a gluon 
scatterer. Such a modification brings with it a change 
in the longitudinal momentum distribution of the gluon 
scatterers which we shall examine in this section. 

We consider the coherent multiple collisions of a gluon 
jet or a gluon jet on n gluons scatterers. From the results 
in Eq. (76), the cross section for the scattering in the 
space-time symmetric and color symmetric state is 



^ ^g^ng 




\&D, (2p • a,)2 dqizdq_^L^^^ 



2a' 



iO 



\Qt\ 



where C^^ is the color factor for the space-time symmet- 
ric and color symmetric component and 



fpi = 4pi • fli 



(85) 



To investigate the longitudinal distribution of the scat- 
terers, we need to investigate the qiz dependence of the 
factor (2p • a^)^ in Eq. (|84| . which can be written as 

(2p-a,)2 = [2p- (2a, + g0]'- (86) 

Because of the A{2p-qi) constraint, the above result gives 

(2p • a,)2A(2p • q,) ^ (4p • a,fA{2p ■ q,), (87) 

where the factor (4p • a^)^ is independent of the longitu- 
dinal recoil qiz of the scatterer. It is convenient to use 



the transfer rapidity to represent the longitudinal mo- 
mentum transfer g^^, 

qiz = mgT sinh^i, (88) 
a'io = \J m^gT + {a^,+ irigT sinh ) ^ , (89) 

where for simplicity the final transverse masses of the 
scatterers are taken to be their average value irigT- To 
make the problem simple, we can take Oiz to have its 
average value (ai^), which is zero in the medium ccnter- 
of-momentum frame. In that approximation, we obtain 
from Eqs. (|HH1) and 



dqiz 

— d^i 

a n 



From Eq. (|84|) . the cross section becomes becomes 
d£,id^2--d(,ndqiTdq2T--dqnT 



where Cg is 



(2aX ffj 16A (2p ■ -ad' 
4 UJ Ui /p. 2 



(90) 



(91) 



(92) 



The quantity Cg is a weak function of S.i and q^rp and can 
be approximated as a constant. The probability distri- 
bution is then a flat function of . The transfer rapidly 
has the upper limit 



^imax cosh 



I Pq 



nnigT 



(93) 



Then as far as is concerned, the average distribution 
is 



dP 1 



d^i £,i: 

The average (qiz) is 



e(e 



po/n - vfigT 

^max 



(94) 



(95) 



With Pq — 10 GeV/c and nigT — 0.6 GeV, we find for 

n = 6, 



-0.6 GeV/c, 
and for for n = 2.4, we find 

q^z 1.4 GeV/c. 



(96) 



(97) 



These estimates indicate that the average longitudinal 
momentum kick acquired by a gluon scatterer is slightly 
smaller than that acquired by a fermion scatterer. They 
are approximately inversely proportional to the number 
of scatterers in a coherent collision. 
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X. SIGNATURES OF BOSE-EINSTEIN 
INTERFERENCE AND COMPARISON WITH 
EXPERIMENTAL DATA 

The results in the above sections provide informa- 
tion on the signatures for the occurrence of the Bose- 
Einstein interference in the coherent collisions of a jet 
with medium partons: 

1. The Bose-Einstein interference is a quantum many- 
body effect. It occurs only in the multiple collisions 
of the fast jet with two or more scatterers. There- 
fore there is a threshold corresponding to the re- 
quirement of two or more scatterers in the multiple 
collisions, n> 2. 

2. Each scatterer has a transverse momentum distri- 
bution of the type l/lgj^j^, which peaks at smah 
values of IqtI- 

3. Each scatterer acquires a longitudinal momentum 
kick that is of order l/2n fraction of the incident 
jet momentum along the incident jet direction . 

4. As a consequence, the final effect is the occurrence 
of collective recoils of the scatterers along the jet 
direction. 

To inquire whether Bose-Einstein interference may cor- 
respond to any observable physical phenomenon, it is 
necessary to identify the scatterers to separate them from 
the incident jet in a measurement. Such a separation 
is indeed possible in Ac/j-Arj angular correlation mea- 
surements of produced pairs with a high-pr trigger [l|]- 
[sot . Particles in the "ridge" part of the correlations with 
IA77I > 0.6 at RHIC and |A?7| > 1.0 at LHC with A0 - 
can be identified as belonging to the medium partons 
because it was observed at RHIC that 

1. The yield of these ridge particles increases approxi- 
mately linearly with the number of participants Q . 

2. The yield of these ridge particles is nearly in- 
dependent of (i) the flavor content, (ii) the me- 
son/hyperon character, and (iii) the transverse mo- 
mentum pT (above 4 GeV) of the jet trigger [1,|3,[1]. 

3. The ridge particles have a temperature (inverse 
slope) that is similar (but slightly higher) than that 
of the inclusive bulk particles, but lower than the 
temperature of the near-side jet fragments 0] 

4. The baryon/mcson ratio of these ridge particles is 
similar to those of the bulk hadrons and is quite 
different from those in the jet fragments [loj . 

With the scatterers as ridge particles that can be sepa- 
rated from the incident high-px jet, the occurrence of the 
Bose-Einstein interference will be signaled by Item (4) of 
the collective recoils of the scatterers (the ridge particles) 
along the jet direction. Such collective recoils will lead 



to the A(j) ^ correlation of the ridge particles with the 
high-pT trigger, as has been observed in angular corre- 
lations of produced hadrons in AuAu collisions at RHIC 
[H-dl, and in pp and PbPb colhsions at LHC 
The collective recoils of the kicked medium partons have 
been encoded into the longitudinal momentum kick (qiz) 
of the momentum kick model that yields the observed 
A(f), A77, and pt dependencies of the angular correlations 
in AuAu collisions at RHIC [3ll - [37j . and pp collisions at 
LHC [13. 

It is of interest to examine Item (3) of the signature 
of the Bose-Einstein interference with regarding to the 
relationship between the (average) magnitude of the lon- 
gitudinal momentum kick, (qiz), and the (average) num- 
ber of scatterers, (n), when such a collective momen- 
tum kick occurs. For the most central AuAu collisions 
at y^SNN = 200 GeV at RHIC, we previously found that 
(/jj)(n) ^ 3.8 where n is the number of kicked medium 
partons and (/j^) is the average attenuation factor for 
the kicked partons to emerge from the collision zone [s^ . 
The value of (ffj) is not determined but a similar at- 
tenuation factor (/j) for jet fragments is of order 0.63 
[33! . We can therefore estimate that for the most central 
AuAu colhsions at ysJvjv = 200 GeV at RHIC, (n) ~ 6. 
For an incident jet of p^^lO GcV/c [2^, the estimates of 
Eq. ([51]) and (72) give 

J0.83GeV/c for a quark scatterer, 
^ \ 0.63 GeV/c for a gluon scatterer. ^ ^ 

These estimates of the momentum kick are of the same 
order as the value of {qtz) ^ 1 GeV/c estimated in [33j 
and 0.8 GeV/c in [3^1, obtained in the momentum kick 
model analysis. 

In another momentum kick model analysis for the high- 
est multiplicity pp collisions at y/s^lv — 7 TcV at the 
LHC, we previously found that {fji){n) ^ 1.5 [l^. The 
value of ridge particle attenuation factor {fn) is not de- 
termined but we can use again the similar attenuation 
factor for jet fragments, (/j) of order 0.63 [11], to esti- 
mate (n) ^ 2.4. For an incident jet of 10 GeV/c, the 
estimates of Eqs. ([51]) and (|95|) give the average scatterer 
longitudinal recoil momentum as 

J 2.1 GeV/c for a quark scatterer, 

I 1.4 GeV/c for a gluon scatterer. ^ ' 

These estimates of the momentum kick arc slightly lower 
than but are of the same order as of the same order as 
the value of (qiz) ^ 2 GeV/c inferred from experimental 
data in the momentum kick model analysis [38| . The ex- 
perimental data give a longitudinal momentum transfer 
that is approximately inverse proportional to the number 
of scatterers. 

With regard to Item (1) for the signature for the oc- 
currence of the Bose-Einstein interference, the presence 
of a threshold implies a sudden increase of the ridge yield 
as a function of the number of kicked medium scatterers, 
n, which increases with the increasing centrality, as rep- 
resented either by an increasing number of participants 
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FIG. 5: (Color online) The ridge yield per high-pT trigger 
as a function of the participant number A'part for nucleus- 
nucleus collisions at -^sjv]v=200 and 62 GeV. The solid cir- 
cular points (for AuAu) and the square points (for CuCu) 
are from the STAR Collaboration HUol]. The curves are the 
momentum kick model results of [341 ] modified to include the 
Bose-Einstein interference threshold effect of n > 2. 



function of iVpart and by uncertainties in the estimates 
of the number of scatterers. Although the experimental 
data appear to be consistent with theory, the large error 
bars and the scarcity of the number of data points in the 
threshold regions preclude a definitive conclusion. 
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FIG. 6: (Color online) (a) Preliminary CMS Collaboration 
data of the ridge yield per trigger in the region of 2 < | At^I < 4 
for pp collisions at 7 TeV, for 2 < Py"*^ < 3 GeV, 1 < p^""" < 
2 GeV, as a function of multiplicities A'^ [lO[l| . (b) The slope 
of the ridge yield, ^(Associated Yield)/dAi' as a function of A''. 



or by increasing multiplicities. The slope of the yield 
as a function of centrality will change from zero below 
threshold to infinite over a small range of centralities at 
threshold. We expect that such a sudden threshold will 
be smoothed out by the fluctuations of the participant 
numbers and multiplicities with respect to the number 
of kicked medium partons. However the change of the 
slope of the ridge yield as a function of centralities will 
remain. Hence, the BE threshold effect will show up as a 
change of the slope of the ridge yield as a function of cen- 
tralities, measured by the number of participants or by 
multiplicities. Equivalently, the ridge yield as a function 
of centralities will appear to have a kink near the thresh- 
old region of centrality. A search for a change of the slope 
or a kink in the ridge yield as a function of centralities 
will allow us to find the BE interference threshold. 

We show in Fig. 5 the experimental ridge yield per 
high-pT trigger as a function of A^part for AuAu and CuCu 
collisions at = 200 and 62 GeV at RHIC H We 
also show in Fig. 5 the theoretical yields obtained in the 
momentum kick model [33j . where the ridge yield at the 
most central collision at A^part = 320 was calibrated as 
n = 6 HI] . With such a calibration, the threshold values 
in A^part at which n = 2 can be located and listed in Table 
I, where the 30% theoretical errors arise from the errors 
in measuring the ridge yield at the most central collision 
at A'part = 320. Theoretical ridge yields from the mo- 
mentum kick analysis in Fig. 10 of Ref. [s^, modified 
to include the Bose-Einstein interference threshold effect 
of n > 2, are shown as the solid curves for AuAu colli- 
sions, and as dashed curves for CuCu collisions in Fig. 
5. The theoretical thresholds in Fig. 5 will be smoothed 
out by the fluctuations of the number of scatterers as a 



In Fig. 6(a), we show the CMS preliminary data on the 
near-side ridge yield per trigger as a function of centrali- 
ties, as measured by multiplicities A', for pp collisions at 
yi=7 TeV [l03|. To search for the BE threshold that 
may show up as a change of slope or a kink of the ridge 
yield as a function of multiplicities, we show the slope of 
the associated ridge yield, (^(Associated Yield) /dN as a 
function of A' in Fig. 6(b). The lines in Fig. 6(a) and 
6(b) join the data points to guide the eyes. The CMS 
preliminary data in Fig. 6(a) and 6(b) indicate a sharp 
change of the slope of the ridge yield in the region around 
A' ~ 50 — 70, which may suggest a threshold for the ridge 
yield at around A' = 50 — 70. 

We show CMS preliminary data [lOOf on the near- 
side ridge yield per trigger for PbPb collision at 
y/sNN— 2.76 GeV/c as a function of the number of 
participants A^part in Fig. 7(a), and the corresponding 
d( Associated Yield)/(iA'part in Fig. 7(b). In Fig. 7, we 
also included the minimum-biased pp data point as an 
open circle for A^part = 2 (for very peripheral PbPb col- 
lisions) , at which the ridge yield is zero [H, [13 ■ The 
integrated near-side associated yield for 4 < p^'^ < 6 
GeV/c and 2 < p^"""" < 4 GeV in the ridge region of 
2 < IA77I < 4 appear to have a kink in the ridge yield as 
a function of A'part at (A'part)~30 in Fig. 7(a). The slope 
increases sharply starting at A'part ^ 20 in Fig. 7(b) and 
reaching a large value at A'part ^ 40. Such a behavior 
may suggest the presence of a ridge yield threshold. The 
location of a possible threshold at (A'part) ^30 for PbPb 
collisions appear to be qualitatively consistent with the 
decrease of (A'part) as a function of increasing collision 
energies, as in(iicated in Table I. It will be of interest 
to examine in future work whether the ridge thresholds 
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as suggested by the CMS pp and PbPb data correspond 
quantitatively to the location of n —2 for the onset of the 
BE interference. 



" PbPb ^s=2.76TeV 
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FIG. 7: Preliminary CMS data of the ridge yield per trigger 
in the ridge region of 2 < IA77I < 4 for PbPb collisions at 
yiT^ =2.76 TeV, for 4 < p*^"*^ < 6 GeV, 2 < < 4 GeV, 

as a function of the number of participants A'^part. [lO0|| . 
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FIG. 8: (Color online) STAR Collaboration data of the 2D 
Gaussian distribution amplitude (ridge yield) for AuAu col- 
lisions at y^sjvjv=200 GeV (solid circles) and 62 GeV (open 
circles), with a low-pT trigger p^ - [l7j . The solid and dashed 
curves are lines joining the data points. The dashed-dot and 
dotted curves represent Glauber linear superposition (GLS) 
estimates. 



Threshold effects for the ridge yield (2D Gaussian 
yield) as a function of iVpart have been observed in an- 
other angular correlation measurements with a low-p^ 
trigger from the STAR Collaboration [l3 - [T7j as shown 
in Fig. 8. We note previously that a fast jet parton pos- 
sesses \ow-pT jet fragments and a minimum-pT^-biased 
\ow-pT trigger can also indicate the passage of a fast par- 
ent jet [38|. As a consequence, ridge particles will also 
be associated with a low-pr trigger. The change of the 
slope of the amplitude of the 2D Gaussian distribution 
(the ridge yield) shown in Fig. 8 indicates the presence of 
a threshold for the ridge yield as a function of centrality. 

TABLE I: Comparison of the locations of the theoretical 
threshold, at which the (average) number of scatterer n is 
equal to 2, with the observed experimental threshold [13 - 
uM , for the sudden increase of the 2D Gaussian distribution 
amplitude and width (ridge component) in AuAu collisions at 
^/sTTiv = 200 and 62 GeV in Fig. 8. 



Collision 


-^SJVJV 


Theoretical 


Experimental 


System 


(GeV) 


Threshold iVpart 


Threshold iVpart 


AuAu 


200 


69±21 


58-86 


AuAu 


62 


146±45 


85-122 


CuCu 


200 


51±15 




CuCu 


62 


103±31 





For AuAu at 200 GeV/c, the experimental ridge yield 
threshold occurs at A^part = 58 — 86 (Fig. 8), which 
can be compared in Table 1 with the theoretical ridge 
yield threshold of A^part — 69 ± 21 as estimated in the 
momentum kick model for n = 2. For AuAu at 62 
GeV/c, the experimental ridge yield threshold occurs at 
Npait = 85 — 122 (Fig. 8), which can be compared with 
the theoretical threshold of A^part = 146 ± 45 as esti- 
mated in the momentum kick model for n — 2 (Table I) . 



These comparisons indicate that experimental data with 
\ow-pT trigger are consistent with the presence of ridge 
thresholds as a function of the number of the medium 
scatterers located at n ^ 2, as in the threshold effect in 
Bosc-Einstein interference. 



XI. CONCLUSIONS AND DISCUSSIONS 

Conventional investigations (TTl - fsll loil - fo^ l on scatter- 
ing of a fast particle with medium scatterers use the po- 
tential model in which the scatterers are represented by 
static potentials and the longitudinal recoils of the scat- 
terers are considered as dependent variables that are ap- 
pendages to the deflected motion of the incident particle. 
A general treatment of coherent collisions necessitates 
the use of the longitudinal recoils of the scatterers as 
independent dynamical variables which are however not 
allowed in the potential model. This leads us to forgo the 
potential model and to turn to the use of Feynman am- 
plitudes for the general treatment of coherent collisions. 

In the Feynman amplitude approach for the coherent 
collision of a fast particle on n scatterers, there are n\ dif- 
ferent orderings in the sequence of collisions along the fast 
particle trajectory at which various virtual bosons are ex- 
changed. By Bose- Einstein symmetry, the total Feynman 
amplitude is the sum of the nl amplitudes for all possi- 
ble interchanges of the virtual bosons. The summation 
of these nl Feynman amplitudes and the accompanying 
interference constitute the Bose-Einstein interference in 
the passage of the fast particle in the dense medium. 

Our interest in examining this problem has been stim- 
ulated by the phenomenological successes of the momen- 
tum kick model in the analysis of the angular correlations 
of hadrons produced in high-energy heavy-ion collisions 
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[3ll - l38| . We seek a theoretical foundation for the origin of 
the fongitudinal momentum kick along the jet direction 
postulated in the model. We explore whether such a lon- 
gitudinal momentum kick may originate from a quantum 
many-body effect arising from the Bose- Einstein interfer- 
ence in the passage of a jet in a dense medium. We take 
note of previous results on the Bose- Einstein interference 
in the emission of real photons and gluons in high-energy 
interactions and in the sum of the ladder and cross-ladder 
loop diagrams in the collision of two particles [83l - [90j . 

We find similarly that in the coherent collisions of an 
energetic fermion with n fermion scatterers at high en- 
ergies in the Abelian theory, the symmetrization of the 
Feynman scattering amplitudes with respect to the inter- 
change of the virtual bosons leads to the Bose-Einstein 
interference, resulting in a sharp distributions at p-qi ^ 0. 
Such coherent collisions are in fact a single collision, ty- 
ing the incident fermion with the n fermion scatterers as 
a single unit. There are then 3(n-|- 1) degrees of freedom, 
subject to the constraints of the conservation of energy 
and momentum. As a consequence, all 3n degrees of free- 
dom of the scatterers can be independently varied. The 
probability distribution in these 3n momentum degrees 
of freedom depends on the Feynman amplitudes and the 
phase-space factors. The Bose-Einstein symmetry con- 
straints oi p ■ Qi ~ limit the transverse momentum 
transfers of the scatterers to small values of qtr- The 
longitudinal momenta of the scatterers get their share 
of longitudinal momenta from the jet, resulting in the 
collective recoils of the scatterers along the jet direction. 

For the coherent collision of an energetic parton with 
parton scatterers in non-Abelian cases, we find that the 
complete Bose-Einstein symmetry in the exchange of vir- 
tual gluons consists not only of space-time exchange sym- 
metry but also color index exchange symmetry. Never- 
theless, there is always a space-time symmetric and color- 
index symmetric component of the Feynman amplitude 
that behaves in the same way as the Feynman ampli- 
tude in the Abelian case, in addition to the occurrence of 
space-time antisymmetric and color-index antisymmetric 
components. For the space-time symmetric and color- 
index symmetric component, the recoiling partons be- 
have in the same way as in collisions in the Abelian case. 
There is thus a finite probability for the parton scatter- 
ers to emerge collectively along the incident trigger jet 
direction, each with a significant fraction of the longi- 
tudinal momentum of the incident jet. The collective 
recoils will lead to the ^ correlation of the ridge 
particles with the high-pT trigger. Such a signature of 
the Bose-Einstein interference may have been observed 
in the A0 correlation of the ridge particles with a 
high-pT trigger in the angular correlation measurements 
of pro duced hadron pairs in AuAu collisions at RHIC 
[H-illl and in pp and PbPb colhsions at LHC [lUll. 
The centrality dependence of the ridge yields i n pp and 
PbPb collisions at LHC with high pT trigger [11 [Tooj 
and in AuAu collisions at RHIC with a low-pr trigger 
[l^ - [l7| gives hints of the presence of a ridge threshold. 



as expected in the quantum many-body effect of Bose- 
Einstein interference. 

Our focus at present is on the recoils of the scatter- 
ers. We have been examining the collision between the 
incident particle and the medium scatterers without ra- 
diation. We find that the scatterer recoils ranges from 
quasi-elastic to substantial fraction of the incident inci- 
dent particle momentum. For a given multiple collision 
processes involving the recoils of the scatterers, the ra- 
diative processes will involve addition external legs in the 
Feynman diagram. They are high order in as. Therefore 
as far as the cross sections are concerned, they gener- 
ally occur with less probabilities. However, as far as jet 
energy loss is concerned, these high-order radiative pro- 
cesses can be important in certain kinematic regions. For 
example, if one restricts oneself to quasi-elastic processes 
with little scatterer longitudinal momentum recoils, then 
radiation can take up more of the jet longitudinal mo- 
mentum, even though they occur with a lower probabil- 
ity than quasi-clastic scattering. One therefore envisages 
that in the investigation of the loss of the incident lon- 
gitudinal momentum, radiative energy loss is important 
in the region of small scatterer longitudinal recoils, but 
the relative importance of the radiative energy loss will 
diminish as one moves toward the regions of greater scat- 
terer longitudinal recoils. 

If one restricts oneself to potential models, then the 
longitudinal momentum and energy losses due to scat- 
terer recoils will be small and radiative energy loss be- 
comes more important than quasi-elastic collision energy 
loss. However, for reasons we examine in Section II and 
discussed earlier in this section, the general treatment 
of coherent collisions necessitates the use of the longitu- 
dinal recoils of the scatterers as independent dynamical 
variables, which are however not allowed in the poten- 
tial model. It is necessary to use the Feynman amplitude 
approach to explore the entire domain of longitudinal 
scatterers recoils. We hope to carry out an analysis of 
radiative processes in conjunction with scatterer recoils 
within the Feynman amplitude approach in the future. 

The present work is based on the high-energy limit 
which allows great simplifications of the algebraic struc- 
tures of Feynman amplitudes. These simplified struc- 
tures bring into clear focus the mechanism of the Bose- 
Einstein interference that changes the nature of the dis- 
tribution function. With the mechanism of the Bose- 
Einstein interference well understood, it may be benefi- 
cial in future work to evaluate the Feynman amplitudes 
without resorting to many of the drastic assumptions and 
simplifications - in order to make quantitative compar- 
ison of theoretical predictions with results of the longi- 
tudinal momentum kick quantities extracted from exper- 
imental data. This is particularly important in regions 
of large longitudinal momentum transfers for which the 
high-energy approximation of having p' not greatly dif- 
ferent from p may not hold. 

Coherence or decoherence in high energy processes oc- 
curs when the scattering amplitude consists of many con- 
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tributions and these contributing amphtudcs interfere. 
Because of its general nature, they come in many differ- 
ent processes in different forms. 

The coherence results wc present here have been ob- 
tained for a single jet and interacting particles as plane 
waves. High-energy collisions has a preference for longi- 
tudinal motion over transverse motion. Therefor, for a 
single jet, as in the case of our main interest, the interfer- 
ence of Feynman diagrams has relevance only with regard 
to longitudinal coherence length along the longitudinal 
direction. For such a configuration, the uncertainties of 
the vertices of boson emission along the longitudinal di- 
rection leads to the coherence in the collision process. 
Uncertainties in the transverse direction come from the 
transverse positions of the scatterers. Such uncertainties 
can be includ ed b y treating the scatterers in terms of 
wave packets 101 1, form factors |l02| . or scatterer wave 
function . In the special case of a purely elastic scat- 
tering of an incident fast particle with medium particles 
(in a nucleus), there can be constructive interference of 
the scattering amplitudes from different particles, leading 
to a sharp diffractivc ela stic peak at the forward angle, 
as described by Glauber jl02l | . 

However, there are circumstances in which interesting 
transverse uncertainties arises, and the amplitudes from 
different transverse sources interfere. An interesting case 
of transverse interference can be found in connection with 
the coherent emission of a gluon from a system of a quark 
jet and an antiquark jet [93. [93j. This case involves the 
propagation of two jets with an opening transverse angle 
9qq and differs from our present case with the propagation 
of a single jet. The amplitude for gluon emission from the 
quark q interfere with the amplitude for gluon emission 
from the antiquark q. As a consequence, in the absence of 
medium, the coherent interference in the transverse direc- 
tion lead to the angular ordering of the emitted gluon in 
such a way that the emitted gluon lies within the opening 
angles of the q-q pair, Ogg . In the presence of a medium, 
as shown in Ref. (92l . l93l | , the quark and the antiquark 
are subject to multiple scatterings along their trajecto- 
ries that change their momenta, their distributions, and 
their propagating phases. These changes of the quark 
and antiquark sources of the gluon due to the medium 
will diminish the strength of the transverse coherent in- 
terference of the emitted gluon that occurs in the absence 
of the medium, leading to a gradual decoherence of the 
gluon emission as a function of increasing medium density 
[9^ [93}. Recent experimental findings on r econ structed 
jets in nuclear collisions at the LHC [l03l Il04| suggest 
that such medium-induced partial decoherence may be 
an important effect. 

Conventional model analysis of collision processes (ttI - 
HH, I9ll - l93 | assumes the ladder- type diagrams such as Fig. 
1(a) and Fig. 2(a), while the cross Feynman diagrams 
such as Fig. 1(b) and Figs. 2(b)-2(f) have been ignored. 
As the derivations in Section HI demonstrate, the con- 
tributions from different cross diagrams destructively in- 
terfere with the ladder Feynman diagram and there is a 



high degrees of cancellation. In general, the neglect of 
these destructive interference by including only the lad- 
der Feynman diagram cannot be justified from a math- 
ematical view point. Only in the special and restrictive 
case of quasi-elastic scattering with the fast particle as- 
sumed to be nearly on the mass shell and with the ad- 
ditional assumption of qio = is it justified to include 
only the ladder diagrams, at the expenses of precluding 
the exploration into other regions of scatterer recoils. 

Many models have bee n prop osed to explain the 
ridge phenomenon (sil - fssl l40l - [7l |. They include the 
collision of jets with medium partons (sil - fssl l60l . [6l| . 
flows and hydrodynamics with initial state fluctuations 

SF 



41 



mm 



|49| . I57l l59l . I6a |. col or-g lass condensate 



58[ , modeling pQCD |62h65| | , parton cascade [6/ 
gluon bremsstrahlung in string formation [68j . strong- 



coupling AdS / CFT [70| , quantum entanglement [72| , and 
BFKL evolution and beyond (T^I • There are however two 
difficulties associated with these phenomenological mod- 
els. Almost all models deal with fragmented parts of the 
data and all models contain implicit and explicit assump- 
tions. 

In the presence of a large number of models and the 
above difficulties, progress can proceed in three fronts. 
First, the models need to cover an extensive set of exper- 
imental differential data over large phase spaces, ccntral- 
ities, and energies, from many different collaborations. 
Secondly, the assumptions of the models need further 
theoretical and observational investigations from funda- 
mental viewpoints. Finally, experimental tests need to 
be proposed to distinguish different models. 

With regard to the first of these three tasks, the mo- 
mentum kick model gives reasonable descriptions for an 
extensive set of differential data of the ridge yield, over 
an extended range of transverse momenta, azimuthal an- 
gles, pseudorapidity angles, centralities, and collision en- 
ergies, from the STAR Collaboration [iHll, the PHENIX 
Collaboration H^-El, the PHOBOS Collaboration [H, 
and the CMS Collaboration [2^. In these momentum 
kick model analyses, a collective longitudinal momentum 
kick on the medium scatterers along the jet direction is 
a central ingredient leading to the successful descriptions 
of the large set of experimental data. Additional analyses 
will be continued to extend the momentum kick model 
to include the effects of collective flows and to cover a 
larger set of new data as they become available. 

With regard to the second of the three tasks, the mo- 
mentum kick model contains the basic assumption that 
the A77 ridge arises mainly from the initial rapidity distri- 
bution of partons prior to the jet collision. Such a basic 
assumption has been examined from the viewpoint of the 
Wigner function of produced particles in a fundamental 
quantum theory of particle production [sl] . Another ba- 
sic assumption concerning the longitudinal nature of the 
momentum kick is now being examined in the present 
manuscript. 

With regard to third of the three tasks, the present 
analysis for jets interacting with medium partons re- 
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veals that the longitudinal momentum kick is a quan- 
tum many-body effect that contains a threshold, requir- 
ing the multiple collisions of the jet with at least two 
partons. This may be in agreement with the onset of 
the ridge yield as a function of the centrality, observed 
by the CMS Collaboration for pp and PbPb collisions 
at LH CjlOOt . and by the STAR Collaboration at RHIC 
0, [T^18| . as discussed in Section X. Furthermore, in 
the momentum kick model with longitudinal momentum 
kicks, the kicked medium partons from back-to-back jets 
possess a {pit,P2t) correlation for both the near side and 
the away side, as observed by the STAR Collaboration 

MM- 

It would be of interest to see how other proposed mod- 
els fare with the three tasks at hand. It will also be 
of interest to see whether they contain the features of a 
rapid rise of the ridge yield as a function of the central- 
ity, and the {pit,P2t) correlation for both the near side 
and the away side. The search for tests to distinguish 
different models will be an on-going research activity. 

Whatever the theoretical descriptions, jets of order 10 
GeV and below are known to be present in high-energy 
collisions pol - l2^ [75| . These jets collide with medium 
partons. They are not hydrodynamical flows, but they 
contribute to the azi muthal a nisotropy and the azimuthal 
Fourier coefficients p5l - [l07t . as well as to two-particle 



correlations. They must be taken into account or be 
subtracted from experimental data in theoretical mod- 
els that do not include the jet effects explicitly. 

In summary, Bose-Einstein interference in the passage 
of a jet in a dense medium is a quantum many-body ef- 
fect that occurs quite generally in a coherent multiple 
collision process, with the threshold of more than two 
scatterers. The manifestation of the Bose-Einstein inter- 
ference effect as collective recoils of the scatterers along 
the jet direction may have been experimentally observed 
in the A(j>^0 correlation of hadrons associated with a 
high-pT trigger in high-energy AuAu collisions at RHIC, 
and pp and PbPb collisions at LHC. The experimental 
observation of ridge thresholds as a function of centrality 
may lend support for the occurrence of the Bose-Einstein 
interference thresholds in the passage of a jet in a dense 
medium. 

Acknowledgment 

The author wishes to thank Profs. C. S. Lam, H. W. 
Crater, Jin-Hee Yoon, Vince Cianciolo, A. V. Koshelkin, 
and Wei Li, for helpful discussions. This research was 
supported in part by the Division of Nuclear Physics, 
U.S. Department of Energy. 



[1] J. Adams et al. for the STAR Collaboration, Phys. Rev. 

Lett. 95, 152301 (2005). 
[2] J. Adams et al. (STAR Collaboration), Phys. Rev. C 

73, 064907 (2006). 
[3] J. Putschke (STAR Collaboration), J. Phys. G34, S679 

(2007). 

[4] J. Bielcikova (STAR Collaboration), J. Phys. G34, S929 
(2007). 

[5] F. Wang (STAR Collaboration), Invited talk at the Xlth 
International Workshop on Correlation and Fluctuation 
in Multiparticle Production, Hangzhou, China, Novem- 
ber 2007, [arXiv:0707.0'8T5] . 

[6] J. Bielcikova (STAR Collaboration), Phys.G34:S929- 
930,2007; J. Bielcikova for the STAR Collaboration, 
Talk presented at 23rd Winter Workshop on Nuclear 
Dynamics, Big Sky, Montana, USA, February 11-18, 
2007, ;arXiv:0707.3100 ; J. Bielcikova for the STAR 
Collaboration, Talk presented at XLIII Rencontres 
de Moriond, QCD and High Energy Interactions, La 
Thuile, March 8-15, 2008, ;arXiv:0806.2261 . 

[7] B. Abelev (STAR Collaboration), Talk presented 
at 23rd Winter Workshop on Nuclear Dynamics, 
Big Sky, Montana, USA, February 11-18, 2007, 
(arXiv:0705.3371 . 

[8] L. Molnar (STAR Collaboration), J. Phys. G 34, S593 
(2007). 

[9] R. S. Longacre (STAR Collaboration), Int. J. Mod. 
Phys. E16, 2149 (2007). 
[10] C. Nattrass (STAR Collaboration), J. Phys. G 35, 
104110 (2008). 

[11] A. Feng, (STAR Collaboration), J. Phys. G 35, 104082 



(2008) . 

[12] P. K. Netrakanti (STAR Collaboration) J. Phys. G 35, 
104010 (2008). 

[13] O. Barannikova (STAR Collaboration), J. Phys. G 35, 
104086 (2008). 

[14] M. Daugherity, (STAR Collaboration), J. Phys. G 35, 
104090 (2008). 

[15] D. Ray, in Talk presented at Tamura Sym- 
posium on Heavy Ion Physics, the Univer- 
sity of Texas at Austin, November 2022, 2008, 
|http:/ /www. ph. utexas.edu/~molly /tamura/ 

[16] D. Kettlcr, (STAR Collaboration), Euro. Phys. Jour. C, 
62, 175 (2009). 

[17] G. Agakishiev et al., (STAR Collaboration), 
larXiv: 1109.4380 (2011). 

[18] T. A. Trainor Phys. Rev C78, 064908 (2008); T. A. 
Trainor and D. T. Kettler, Phys. Rev. D 74, 034012 
(2006); T. A. Trainor, Phys. Rev. C 80, 044901 

(2009) ;T. A. Trainor, J. Phys. G 37, 085004 (2010). 
[19] M. van Leeuwen, (STAR Collaboration), Eur. Phys. J. 

C 61, 569 (2009). 
[20] A. Adare, et al. (PHENIX Collaboration), Phys. Rev. 

C 78, 014901 (2008). 
[21] M. P. McCumber (PHENIX Collaboration), J. Phys. G 

35, 104081 (2008). 
[22] Chin-Hao Chen (PHENIX CoUaboration), "Studying 

the Medium Response by Two Particle Correlations", 

Hard Probes 2008 Intern. Conf. on Hard Probes of High 

Energy Nuclear Collisions, A Toxa, Galicia, Spain, June 

8-14, 2008. 

[23] Jiangyong Jia, (PHENIX CoUaboration), J. Phys. G 35, 



21 



104033 (2008). 
[24] M.J. Tannenbaum, Eur. Phys. J. C 61, 747 (2009). 
[25] E. Wenger (PHOBOS Collaboration), J. Phys. G 35, 

104080 (2008). 

[26] CMS Collaboration, JHEP 1009, 091 
(2010), arXiv:1009.4122 ^ 

[27] CMS Collaboration, ar Xiv:1105.2"438l (2011). 

[28] CMS Collaboration, arXiv:l20i:3158 (2012). 

[29] ATLAS Collaboration, arXiv: 1203.3087 (2012). 

[30] ATLAS Collaboration, Phys. Lett. B708, 249 (2012). 

[31] C. Y. Wong, Phys. Rev. C 76, 054908 (2007). 

[32] C. Y. Wong, Chin. Phys. Lett. 25, 3936 (2008). 

[33] C. Y. Wong, J. Phys. G 35, 104085 (2008). 

[34] C. Y. Wong, Phys. Rev. C 78, 064905 (2008). 

[35] C. Y. Wong, Phys. Rev. C 80, 034908 (2009). 

[36] C. Y. Wong, Phys. Rev. C 80, 054917 (2009). 

[37] C. Y. Wong, Nonlin. Phenom. Complex Syst. 12, 315 
(2009), arXiv:0911.3583 . 

[38] C. Y. Wong, Phys. Rev. C 84, 024901 (2011). 

[39] C. Y. Wong, Invited talk presented at the 35th Sym- 
posium on Nuclear Physics, Cocoyoc, Mexico, January 
3, 2012, to be published in lOP Conference Series, 
[arXiv: 1203.4441 (2012)] . 

[40] E. Shuryak, F'hys. Rec. C 76, 047901 (2007). 

[41] S. A. Voloshin, Nucl. Phys. A 749, 287 (2005). 

[42] C. B. Chiu and R. C. Hwa Phys. Rev. C 79, 034901 
(2009). 

[43] R. C. Hwa and C. B. Yang, Phys.Rev. C 67 034902 
(2003); R. C. Hwa and Z. G. Tan, Phys. Rev. C 
72, 057902 (2005); R. C. Hwa and C. B. Yang, 



[nucl-th/0602024 



[44] C. B. Chiu and R. C. Hwa Phys. Rev. C 72, 034903 
(2005). 

[45] R. C. Hwa, Phys. Lett. B666, 228 (2008). 
[46] V. S. Pantuev, ;arXiv:0710.1882 . 

[47] A. Dumitru, F. Gelis, L. McLerran, and R. Venu- 

gopalan, Nucl. Phys. A810, 91 (2008). 
[48] S. Gavin, and G. MoscheUi, J. Phys. G 35, 104084 

(2008) . 

[49] S. Gavin, L. McLerran, and G. MoscheUi, Phys. Rev. C 
79, 051902 (2009). 

[50] N. Armesto, C. A. Salgado, and U. A. Wiedemann, 
Phys. Rev. Lett. 93, 242301 (2004). 

[51] P. Romatschke, Phys. Rev. C 75 014901 (2007). 

[52] A. Majumder, B. MuUer, and S. A. Bass, Phys. Rev. 
Lett. 99, 042301 (2007). 

[53] A. Dumitru, Y. Nara, B. Schenke, and M. Strickland, 
Phys. Rev. C 78, 024909 (2008); B. Schenke, A. Du- 
mitru, Y. Nara, and M. Strickland, J. Phys. G 35, 
104109 (2008). 

[54] R. Mizukawa, T. Hirano, M. Isse, Y. Nara, and A. 

Ohnishi, J. Phys. G 35, 104083 (2008). 
[55] Jianyong Jia and R.. Lacey, Phys. Rev. C 79, 011901 

(2009) . 

[56] Jianyong Jia, Eur. Phys. J. C 61, 255 (2009). 

[57] Y. Hama, R. P. G. Andrade, F. Grassi, W.- 
L. Qian, Talk presented at ISMD2010, 21-25 
September, 2010, University of Antwerp (Belgium), 
farXiv: 1012. 1342 ; R .P.G.Andrade, F. Gardim, F. Grassi, 
Y.Hama, W.L.Qian [arXiv:1107.0"216i . 

[58] A. Dumitru, K. Dusling, F. Gelis, J. Jalilian-Marian, 
T. Lappi, and R. Venugopalan, Phys. Lett. B697 21 
(2011), larXiv: 1009.5295 . 

[59] K. Werner, lu. Karpenko, K. Mikhailov, and T. 



and Y. 



V. 



de Deus, and 



Kovchegov, 
C. Pajares, 



Pierog "arXiv: 1104.3269] : Fu-Ming Liu, K. Werner, 

arXiv: 1106.5909]. 
[60] R. C. H wa, C. B. Yang, Phys. Rev. C 83, 024911 (2011), 

irXi^:1011.0965 . 
[61] C. B. Chiu and R. C. Hwa, [axriv: 1012:3486]. 
[62] T. A. Trainor, ^arXiv:1008 .4757i T. A. Tr ainor, 

larXiv:1011.635n T. A. Trainor, larXiv:1012.2373l 
[63] T. A. Trainor and D. T. Kettler, Phys. Rev. C 83, 

034903 (2011). 
[64] T. A. Trainor and D. T. Kettler, arXiv: 1010.3048, 
[65] T. A. Trainor and R . L. Ray, ..arXiv: 1105.542.8 ; R. L. 

Rav. larXiv:1106.5023l 

[66] B. Sc henke, [arXiv:1106.6012] : B. Schenke, S. Jeon, C. 

Gale, llarXiv:1109.6289| . 
[67] H. Petersen, C. Greiner, V. Bhattacharya, S. A. Bass, 

[arXiv:1 105.0340 . 
[68] B. A. Arbuzov, E. E. Boos, and V. I. Savrin, 

|arXiv:1104.1283| . 
[69] M. Yu. Azarkin, L M. Dremin, and A. V. Leonidov, 

|arXiv:1102 .3258 . 
[70] IL R. Grigoryan, 

|arXiv :1012.5431 . 
[71] L Bautista, J. Dias 

[arXiv:1011.T870] . 

[72] I. O. Cherednikov and N. G. Stefanis, 'arXiv: 1010.4463] . 
[73] I. M. Dremin, and V. T. Kim, arXiv:'l010.0918]. 
[74] E. Levin, A. H. Rezaeian, arXiv:1105.3275 . 
[75] X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 

(1991); X. N. Wang and M. Gyulassy, Phys. Rev. Lett. 

68, 1480 (1992). 
[76] J. D. Bjorken, Fermilab-Pub-82-059-THY (1982). 
[77] M. Gyulassy and X. N. Wang, Nucl. Phys. B 420, 583 

(1994). 

[78] R. Baier, Yu. L. Dokshitzer, A. H. S. Peigne, and D. 

Schiff, Nucl. Phys. B478, 577 (1996). 
[79] U. A. Wiedemann, Nucl. Phys. B582, 409 (2000). 
[80] M. Gyulassy, P. Levai, and I. Vitev, Nucl. Phys. B594, 

371 (2001). 

[81] M. Djordjevic and M. Gyulassy, Nucl. Phys. A733, 265 
(2004). 

[82] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 
C69, 034910 (2004); J. Adam et al. (STAR Collabora- 
tion), Phys. Rev. Lett. 91, 172302 (2003); B. B. Back 
et al. (PHOBOS Collaboration), Phys. Lett. B578, 297 
(2004); L Arsene et al. (BRAHMS Collaboration), Phys. 
Rev. Lett. 91, 072305 (2003); B. I. Abelev et al. (STAR 
Collaboration), Phys. Rev. Lett. 98, 192301 (2007); S. 
S. Adare et al. (PHENIX Collaboration), Phys. Rev. 
Lett. 97, 252002 (2006). 
[83] H. Cheng and T. T. Wu, Phys. Rev. 186, 1611 (1969). 
[84] H. Cheng and T. T. Wu, Expanding Protons: Scattering 

at High Energies , M. I. T. Press, 1987. 
[85] Y. J. Feng, O. Hamidi-Ravari, and C. S. Lam, Phys. 

Rev. D 54, 3114 (1996). 
[86] Y. J. Feng, O. Hamidi-Ravari, and C. S. Lam, Phys. 

Rev. D 55, 4016 (1997). 
[87] C. S. Lam, Lectures given at the First Asia Pacific 
Workshop on Strong Interactions, Taipei. August 1st 
to 27th, 1996, hep-ph/9704240 
[88] C. S. Lam and K.F. Liu, Nucl. Phys. B 483, 514 (1997). 
[89] C. S. Lam and K.F. Liu, Phy. Rev. Lett. 79, 597 (1997). 
C. S. Lam, Chin. Jour. Phys. 35, 758 (1997), 



[90] 
[91] 



hep-ph/9805210 

G. Ovanesyan and I. 



Vitev , JHEP 1106 (2011) 080. 



22 



[92] Y. Mehtar-Tani, C.A. Salgado, and K. Tywoniuk, 

J.Phys. G38, 124063 (2011). 
[93] N. Armesto, Hao Ma, Y. Mehtar-Tani, C. A. Salgado, 

K. Tywoniu, JHEP 1201, 109 (2012), 
[94] R. J. Glauber, in Lectures in Theoretical Physics, edited 

by W. E. Brittin and L. G. Dunham (Interscience, N.Y., 

1959), Vol. 1, p. 315. 
[95] In this manuscript, we use the Feynman rule in Ref. [13] 

which includes explicitly an overall multiplicative phase 

factor of (— i). Such a overall phase (— i) is not explicitly 

stated but implicitly implied in Appendix A-4 of Ref. 

[96] C. Itzykson and J.B. Zuber, Quantum Field Theory, Mc- 

Graw HiU, (1980). 
[97] V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, 

Quantum Electrodynamics, Perganion Press, 1982. 
[98] C. Y. Lo and H. Cheng, Phys. Rev. D13, 1131 (1976); 

P. S. Yeung, Phys. Rev. D 13, 23062317 (1976). 
[99] H. Cheng, J. A. Dickinson, and K. Olaussen, Phys. Rev. 



D 23, 534 (1981). 

[100] CMS CoUboration, in CMS Physics Anal- 
ysis Summary HIN- 11-006, available on 

the CERN CDS Information Server at 

http://cdsweb.cern.ch/record/1353583/files/HIN-ll-006-pas.pdf 

[101] M. E. Peskin and D. V. Schroeder An Introduction 
To Quantum Field Theory, Addison- Wesley Publishing 
Company, (1995). 

[102] R. J. Glauber, in High-Energy Physics and Nuclear 
Structure, edited by S. Devons, Press Press, N.Y. , 1970. 

[103] G. Aad et al. (Atlas Collaboration), Phys. Rev. Lett. 
105, 252303 (20101. [arXivilOll. 61821 

[104] S. Chatrchyan et al. (CMS Collaboration), Phys. Rev. 
C84 , 024906 (2011). [arXiv:1102. 19571 

[105] C. Y. Wong, Phys. Lett. B 88, 39 (1979). 

[106] J. Y. OUitrauh, Phys. Rev. D 46, 229 (1992). 

[107] S. Voloshin, Y. Zhang, Z. Phys. Rev. C 70, 665 (1996). 



